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ABSTRACT 


A mathematical model and computer code have been developed 
to compute the reliability of a monopropeliant blowdown spacecraft orbit 
adjustment propulsion system as a function of time. The reliability model 
interfaces with a computer code that models the performance of a blowdown 
(unregulated) monopropeliant auxiliary propulsion system. The computer 
code acts as a performance model and as such gives an accurate time his- 
tory of the system operating parameters. The basic timing and status 
information is passed on to and utilized by the reliability model which 
establishes the probability of successfully accomplishing each required 
orbit adjustment. 

The mathematical model and computer code were developed as a 
background effort to verify the concepts prior to writing a reliability 
model based on the current NERVA propulsion system. 
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RELIABILITY MODEL 
OF A 

MONOPROPELLANT AUXILIARY 
PROPULSION SYSTEM 

by Joel S. Greenberg 


INTRODUCTION 

A measure of how well a system performs or meets its design 
objectives is provided by the concept of system reliability. In general, 
reliability can be defined as the probability of successful system oper- 
ation in the manner and under the conditions of intended use. Since space 
systems are normally designed to achieve multiple time dependent objec- 
tives, it is important to establish the reliability of achieving these 
objectives so that appropriate design trade-offs can be made. In order 
to demonstrate the importance of evaluating the reliability of multiple 
time dependent objectives and the basic techniques employed, a relatively 
simple mission was considered. Specifically, a mathematical model and 
associated computer code has been developed which computes the reliability 
of a monopropellant blowdown hydrozine spacecraft auxiliary propulsion 
system as a function of time. The propulsion system is used to adjust or 
modify the spacecraft orbit over an extended period of time. The multiple 
orbit corrections are the multiple objectives which the auxiliary propul- 
sion system is designed to achieve. Thus the reliability mode] computes 
the probability of successfully accomplishing each of the desired orbit 
corrections. To accomplish this, the reliability model interfaces with a 
computer code that models the performance of a blowdown (unregulated) mono- 
propellant auxiliary propulsion system. The computer code acts as a per- 
formance model and as such gives an accurate time history of the system 
operating parameters. The basic timing and status information is passed 
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on to and utilized by the reliability model which establishes the proba- 
bility of successfully accomplishing the orbit corrections. 


II. GENERAL PROPULSION SYSTEM DESCRIPTION 

Monopropellant auxiliary propulsion systems are generally low 
thrust systems which are used for spacecraft attitude control and station 
keeping. Figure 1 illustrates the blowdown hydrazine propulsion system 
which has been modeled. The system consists of two identical monopropel- 
lant half-systems which are interconnected through a valve system such 
that the propellant from either half -system can be expended via either 
thruster. The half-system interconnecting valve system is normally in 
the closed position. The initial pressure in the propellant tank is such 
that the ullage gas will expand (blow-down) against the propellant, forcing 
it through the system during thrust periods until the propellant is depleted. 
The blow-down mode may be staged. That is, a pressurant tank may be placed 
in series with the main propellant tanks and at a predetermined propellant 
tank pressure the associated squib valve opens and the system is repres- 
surized (i.e., blowdown occurs). Thus, sufficient pressure is maintained 
during normal operation to ensure that the desired amount of propellant can 
be expended. 

The general propulsion system configuration which is modeled is as 
indicated in Figure 1. Within this configuration a number of alternatives 
are possible. These alternatives may be summarized as follows: 

(a) A single stage blow-down system may be considered. This 
implies the elimination of the high pressure system and 
its associated squib valve. 
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(b) All squib valves are considered to consist of n valves in 
parallel . 

(c) The high pressure system can consist of m high pressure tanks 
in a parallel arrangement. 

(d) The half-system interconnecting valve system can consist of 
either squib or solenoid valves . 

(e) Either a single fill valve can be used for each pressurant 
and propellant tank or a single propellant fill valve can be 
used for each half-system. 

(f) Solenoid valve systems can consist of either a single valve, 
dual series valves, dual parallel valves, quad or quad con- 
nected valves . 

Ill . PERFORMANCE MODEL* 

The performance model and associated computer code (the entire com- 
puter code for both the performance and reliability models is contained in 
Appendix I) yields an accurate time history of the propulsion system operating 
parameters. The mathematical equations used to represent the components are 
the basic orifice flow equation, thermodynamic expansion equation, and the 
various equations relating the rocket engine parameters . In order to run 
the performance program, the system hardware must be completely defined. 

The basic performance input data is summarized in Table I. Additional in- 
formation describing thruster characteristic performance curves and specific 


*The performance model was developed by Ralph Lake of RCA's Astro Electronic 
Division as part of course (E581- 1970/71) requirements at Princeton Univer- 
sity. 
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TABLE I 


Input Variable 

VEJIN 

Z 

CFIN 

DEN 

AT 

PFIL 

STG 

CD 

A 

P 

DPT 

PMO 

VTK 

VTINT 

PW 

WSC 

L 

ULV 

DELV 

PRES 

PCMIN 

PCIN 

DP 


PERFORMANCE INPUT DATA 

Description 

Effective jet velocity (m/sec) 

Compressibility factor 

Nozzle thrust coefficient , 

3 

Propellant density (kg/m ) 

2 

Nozzle thrust area (m ) 

2 

Initial tank pressure (n/m ) 

Number of blowdown stages 

Orifice coefficient (all components) 

Orifice area (all components) 

Gas expansion exponent 

2 

>. Initial tank internal pressure drop (n/m ) 

Initial propellant mass (kg) 

3 

Tank volume (m ) 

3 

Volume of internal tank construction (m ) 

Pulse width (sec) 

Spacecraft mass (kg) 

Number of maneuvers to be considered 

3 

.... High pressure ullage volume (m ) 

Velocity increment for each maneuver (m/sec) 

2 

Repressurization point pressure (n/m ) 

2 

Minimum operating chamber pressure (n/m ) 

2 

Base point chamber pressure (n/m ) 

System tubing diameter (m) 
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component parameters must also be provided (for example, effective jet 
velocity vs. chamber pressure, nozzle thrust coefficient vs. chamber pres- 
sure) . 

The performance model deals with the flow restrictors (valves, ori- 
fice, filters) and sizes the trim orifice to achieve the desired base point 
parameters for beginning of life operation. For a given mission, the velocity 
increment schedule, the thruster pulse width and number of blowdown stages 
are the primary inputs. Efficiency factors for the thruster depend on both 
operating pressure and the pulse width. These factors are introduced during 
an iteration process in terms of the effective jet velocity versus chamber 
pressure and thruster efficiency versus pulse number.* The model iterates 
the pulsing performance until the required velocity increment for that (Ktfe) 
maneuver is reached. The model keeps track of the number of impulses re- 
quired to accomplish the manuever. 

i 

The model iterates the gas expansion process with time (in terms of 
the blowdown conditions) to determine the propellant flow rate and calculates 
tank pressure, chamber pressure, propellant remaining, effective jet velocity, 
specific impulse, thrust, mass flow rate, total impulse, and burn time. , When 
the propellant mass decreases to zero the iteration terminates. The program 
may also be terminated if the system goes below the minimum operating cham- 
ber pressure. All of the above parameters are printed out for each of the 
maneuvers . 


*A11 curves are input in table form and a subroutine (STINT) is utilized to 
interpolate and return the correct values to the main program. 


The output (see Appendix II) lists the important system operating 
parameters as a function of the maneuver. This allows the propulsion system 
performance to be monitored as a function of time and indicates whether or 
not the modeled system fulfills the mission requirements. 

t 

RELIABILITY MODEL < 2-7 ) 

General 

The term reliability, as used herein, denotes the probability that a 
component, subsystem, or system (as the case may be) will perform its intended 
functions adequately under defined operating conditions at a designated time 
for a specified operating period. Reliability predicts mathematically the 
equipments behavior under expected operating conditions. Reliability is syn- 
onymous with probability of survival or probability of success. 

The starting point of reliability analysis is the determination or 
specification of the intended functions and successful operation or adequate 
performance (i.e., the definition of success) . To judge adequate performance 
involves the observation of inadequate performance in operation', therefore, 
the observation of malfunctions and failures which violate the requirement 
for an "adequate performance." The frequency of failures and malfunctions 
is an important parameter used in the mathematical formulation of reliability. 
This parameter is referred to as the failure rate . It is usually measured in 
terms of the number of failures per unit time or cycle of operation. 

The reliability model considers both random or chance failures and wear- 
out failures. Random or chance failures are characterized by a failure rate 
which is independent of time (i.e., a constant). Under this condition the 



probability distribution of time to failure is given by 

f (t) = A exp (- A t) 

where the failure rate A may be statistically estimated by experiment. 
Reliability , that is the probability that the device is functioning properly 
at time t, is given by 

R (t) = 1 ~ / f(t) dt = exp (-At) 

where A is the constant chance failure rate and t is the operating time 
for which it is desired to determine the reliability, R(t). The exponential 
distribution is independent of device age. When the number of operating 
cycles is more meaningful than time,- the reliability equation can be written 
as 


R(c) = exp (-At) 

where R(c) is the reliability at c operating cycles, or the probability of 
survival through c cycles. 

The normal density function is used to characterize wearout failures 
and is of the form 

(T - M) 2 1 


where T is the component or device age, M is the expected or mean life, and 

a is the standard deviation of the lifetimes. The reliability or probability 

of survival to time T is thus given by 

t 00 2 

R(T) = 1 - ft (t) dt = 1 f exp[ - ^ ] dt 

y/2* <* 2 <r 2 

Since chance and wearout phenomena usually occur simultaneously, it is 
necessary to evaluate their combined; effects . Thus it is desired to evaluate 


f(T) = 1 exp 


9 



the reliability of a mission having a duration of t hours and employing 
a component that has an age of T hours at the beginning of the mission. 

The combined probability of failure in t equals the probability of failing 
of chance and /or wearout in the interval t, at the age of T. Therefore, 

Q(t) - Q c (t) + F w (t) - Q c (t) • F w (t) 

where Q(t) is the probability of failing at time t, Q c (t) is the probability 

of failing at time t due to chance failures, and F w (t) is the probability of 

failing at time t due to wearout. F w (t) is the a posteriori or conditional 

probability of a failure given survival to an age T. Therefore, it can be 
( 2 ) 

shown v ' that the probability of successful operation is given by 


_ R(t) 


- Xt 


-(T-M) 2 /2a 2 

T^t 6 dT 


T 


/V (T - M) Z /2o ‘ dT 


A similar expression can be written in terms of cycles of operation. 

System reliability is normally a calculated or computed quantity. It 
is based upon consideration of the components used in the system, how they 
are used, their modes of failure, and their probability of successful opera- 
tion. Component reliabilities are normally obtained from tests which yield 
information about failure rates. System reliability is thus an extrapolation 
of the component reliabilities using a mathematical model to describe the 
system operation. 

From basic probability theory it can be seen that the probability of 
successful operation of n components connected in series is 

n 

R g (t) = R 1 (t)-R 2 (t)***-R n (t) = exp (-ZA.t) 

* ~ " i=l 

This is called the "product law of reliabilities" for components operating in a 
serial arrangement. A similar equation can be developed for components operating 
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in parallel (i.e., in order to have a system failure all parallel components 
must fail). Therefore the probability of no failures is 

Qp(t)=Q 1 (t)‘Q 2 (t)>-"Q n (t)= i n i Q.(t) 
where Q^(t) = probability of no failure = 1 - R^(t). This is termed the "pro- 
duct law of unreliabilities in parallel operation." This can also be written 
as 

E p (t) = 1 - 1= n. t 1 - 


where R^(t) represents the probability of successful operation (i.e., at least 
one of the components operating satisfactorily) . 

These are thus the basic reliability concepts which are utilized in the 
reliability model described in the following section. 


Mathematical Model 

The mathematical model to be described in detail computes the reliability 
of the monopropellant blowdown spacecraft auxiliary propulsion system illustra- 
ted in Figure 1. The model establishes the probability of performing each of a 
series of orbit corrections. The model is based upon the following assumptions: 

a. The propulsion system consists of two symmetric half-systems which 
are interconnected through a valve system (normally closed) . The 
interconnecting valve system is opened when one of the components 

in the thruster system (fill valve, pressure transducer , lines , filter, 
orifice, valve system, thruster) in use fails. 

b. The valve systems which control the propellant flow (in series with 

■V 

the propellant tanks) are normally closed and must therefore open 
to allow propellant flow. 

c. Each half-system contains three propellant tanks with associated 
valve systems. All three valves are normally closed and are opened 
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simultaneously to allow propellant flow. Therefore, when there is 
no failure, propellant is drawn in equal amounts from the three pro- 
pellant tanks. 

d. The useable propellant of a single half-system is expended before 
the second half— system is utilized. This implies that, repressuriza— 
tion of a single half-system is performed prior to using the second 
half-system. 

e. If a leak develops in the pressurant system, the associated half- 
system propellant system fails. 

Definitions of all input variables (for the reliability model) are given 
in Table II, 

As discussed previously, the reliability model is designed to interface 
with the performance model. This imposes several constraints upon the reliabil- 
ity model. The performance model is implemented such that the performance of 
each half-system is evaluated independently. Thus the performance model must 
be cycled through twice - once for each haif-system. When the first half-system 
is being evaluated HALFZ = 1 and when the second half-system is being evaluated 
HALFZ - 2. 

Since the total propellant mass is distributed equally between the two 
half-sys terns , 

PMO = PMO/ 2 

When HALFZ =1 

Then WSC = as input 

PM1 k = PM k 

PM2, = PMO 
k 
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TABLE II 


RELIABILITY MODEL - INPUT DATA 


Term Dimension 


Max. Value 


Definition 


TIMEZ (K) 100 


XXXX.XX Time (days) of occurance of K orbit 
correction. 


NBZ 


1 


XX 


Average number of brazed joints per 
thruster engine assembly. 


NPZ 


1 


XX Two symetric half-systems are considered 

with each half system consisting of 
three propellant tank groups. NPZ rep- 
resents the number of propellant tanks 
per group . 


NS1Z 1 


XX 


Number of parallel squib valves in 
high pressure system. 


NS2Z 1 


XX 


NSCZ 1 XX 


Number of parallel squib valves in 
half-system connecting value system . 
When CONVZ=SOL then set NS2Z=0. 

r 

Number of catalyst screens per thruster. 


NTZ 1 XX Number of high pressure tanks per 

half -system. 

NWHPZ 1 XXX Number of weld connections (per half- 

system) between the high pressure 
squib valve and the pressurant and 
propellant tanks. 

NWPZ 1 XXX Number of weld connections (per prop- 

ellant tank sys.) between propellant 
tank and propellant tank solenoid 
valve system. 

NWTZ 1 XXX 


Number of weld connections (per half- 
system) between propellant tank solenoid 
valve systems and thruster. 
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TABLE II 


Term 

VSC 

VSP 

VST 


CONVZ 

SPFVZ 


MBZ 

STDZ 

PFS£ 

LCZ 


RELIABILITY MODEL - INPUT DATA (contd) 


Dimension Max. Value Definition 

1 it^n* "DP"* Specification of type of valve system 

"Q"» "Q c " configuration to be used for the half- 

system connecting solenoid valve system. 
Only required when CONVZ="SOL". 


DS”, "DP", Specification of type valve system 
QC configuration to be used for the prop- 

ellant solenoid valve system. 


"S'7 "DS", "DP", 


"ntt 


'Q", "QC' 


Specification of type of valve system 
configuration to be used for the 
thruster solenoid value system. 


'SQ" or "SOL" 


Specification of the type of valve 
to be used for connecting the two 
half-systems. "SQ" refers to squib 
and "SOL" refers to solenoid. 


"Y" or "N" 


When SPFVZ = Y, a single propellant fill 
valve will be used for each half-system. 
When SPFVZ=N, a single propellant fill 
valve will be used for each propellant 
tank. 


XXX 


Average or mean number of cycles for 
bladder failure rate due to wear out 
phenomena . 


XXX 


Standard deviation (cycles) of bladder 
wear out failures. 


.XXXXXX 


Probability of a squib firing 
fully when required. 


success- 


1 XXX.XXXX 


Random or chance failure rate of valve 
caps (failures/10 6 hours) . 
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TABLE II 


Term 

LCLZ 

LOPZ 

LVSZ 

LVZ 

LWZ 

LPZ 

LFZ 

LBRZ 

LSCZ 

LBZ 


RELIABILITY MODEL - INPUT DATA (contd ) 

Dimension - Max. Value Definition 

1 XXX.XXXX Randon or chance failure rate of solenoid 

valve for failing in a closed position 
(failure/10^ cycles). 


XXX.XXXX Random or chance failure rate of solenoid 
valve for failing in an open position 
(failures/lCr cycles). 


XXX.XXXX Random or chance failure rate of solenoid 
valve for excessive leakage past the valve 
seat (failures/10^ hours) . 


XXX.XXXX Random or chance failure rate of fill 

valves (quick disconnect) (failures /10^ 
hours) . 


XXX.XXXX Random or chance failure rate of welded 
connections (failures/10^ hours). 


XXX.XXXX 


Random or chance failure fate of pressure 
transducer (failure/10 hours) . 


XXX.XXXX Random or chance failure rate of filter 
media (failure/10^ hours) . 


XXX.XXXX Random or chance failure rate of brazed 
joints (failures/10^ hours). 


XXX.XXXX Random or chance failure rate of catalyst 
envelope screens (failures /106 hours) . 


1 XXX.XXXX Random or chance failure rate of bladder 

(failures/106 hours) . 
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Where PMl^ = PMO 

WSC is the mass of the spacecraft, PM^ is the propellant remaining at the 
completion of the Rt* 1 orbit correction in the half-system under consideration, 
and PMl k and PM2 k represent the propellant in the first and second half- 
systems, respectively. It should be noted that when HALFZ = 2, the space- 
craft mass must be initialized to take into account propellant remaining 
in the first half-system. 

A schedule of velocity increments for each orbit correction is a neces- 
sary input for the performance computations. Similarly, the time (TIMEZ^) of 
each of the K orbit corrections is a necessary input for the reliability com- 
putations. Since the input data is specified in days and the computations are 
performed in terms of hours, 

TIMEZ k = 24.0 • TIMEZ k 

The performance model keeps track of which half-system is in use. The 
determination of the half-system in use is based upon the assumption that a 
half-sys tern will be used until all usable propellant is expended. This in- 
cludes the repressurization cycle. Therefore, 

When first half-sys tem is in use 
Then HALFZ =1 

When second half -system is in use 
Then HALFZ = 2 

It is necessary, as will be seen in the following pages, to have know- 
ledge as to when repressurization takes place. To achieve this the variable 


When HALFZ = 2 

Then WSC = WSC + PM1, 


p M! k = PMl^ 


Prior to output, set PMj^ = PM]^ + PM2 k 
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MODE is defined. MODE is equal to zero prior to repressurization and equal 
to 2 after repressurization. M0DE1 and M0DE2 refer to the repressurization 
status of half-systems 1 and 2 respectively. The value of MODE is available 
from the performance portion of the model. Thus, 

When HALFZ = 1 and MODE = 0 

Then M0DE1 = 0 and MODE 2 = 0 
When HALFZ = 1 and MODE =2 

Then M0DE1 = 2 and M0DE2 - 0 
When HALFZ = 2 and MODE = 0 

Then M0DE1 = 2 and MODE 2 = 0 
When HALFZ =2 and MODE =2 

Then M0DE1 = 2 and MODE 2 = 2 

It is assumed that all fill valves used in the system are capped and 
are of the quick disconnect type. In order to have a failure both the fill 
valve and the cap must fail. Therefore, the reliability of the fill valves, 
RFVZ^, at the completion of the orbit correction is 

RFVZ k - 1.0 - [l.J0 - exp (-LVZ • TIMEZ k • 10~ 6 )] 

• [ 1.0 - exp (-LCZ • TIMEZ k • 10 -6 ] 

where LVZ and LCZ are the chance failure fates of fill valves and valve 
caps (failures per 10 ^ hours) , respectively. 

The reliability of the high pressure tank system, RTZ k , is given by 
RTZ k = exp [ -NTZ * LWZ * TIMEZ k * 10~ 6 ] 

where it is assumed that the failure mode is due to leakage through the 
weld joints. NTZ and LWZ are the number of parallel high pressure tanks 
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per half— system and the random or chance failure rate of welded connec- 
tions (failures per l£f^ hours), respectively. Similarly the reliability of 
the propellant tank weld ring is 

RPTZ fc = exp [ -LWZ TIMEZ fc • 1 gf 6 ] 

The reliability of the squib valve system (SV1) attached to the high 
pressure tanks is given by 

RSV1Z = 1.10 - ( 1.0 - PFSZ) NS1Z 

where PFSZ is the probability of a squib firing successfully when required and 
NS1Z is the number of squibs in parallel. 

There are a number of connections in the system which consist of weld 
(or brazed) joints. The welded connections have been broken down into three 
groups namely (a) those between the squib valve (SV1) and the pressurant and 
propellant tanks, (b) those between the propellant tank and the associated 
solenoid valve system, and (c) those between the propellant tank solenoid 
valve system and the thruster. The reliability of these three groups of 
connections are 

RLHPZ k = exp [ -NWHPZ • LWZ • TIMEZ k • 10~ 6 ] 

RLPPZ k = exp [ -NWPZ • LWZ ' TIMEZ k ' 10~ 6 ] 

RLTZ fc = exp [ -NWTZ * LWZ • TIMEZ fc • 10~ 6 ] 
where NWHPZ, NWPZ, NWTZ are the number of connections i n each of the above 
three groups, respectively. 

The reliability of the half-system connecting valve system depends upon 
whether squib or solenoid valves are used. When squib valves are used, i.e.. 
When CONVZ = "SQ" (or not "SOL") 

Then 

RSV2Z k = 1.0 - ( 1.0 - EFZ) NS2Z 
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where NS2Z is the number of squib valves in parallel and CONVZ is the input 
variable which specifies whether a squib (SQ) or a solenoid (SOL) system is 
to be used. 

The reliability of solenoid valves depends upon both time and the 
number of cycles of operation. The solenoid valve system used for the half- 
system valve requires only one cycle of operation, i.e., it is normally 
closed and is required to open only when the thruster system in use fails. 

The failure modes of a solenoid valve are failure to open when required, 
failure to close when required, and excessive leakage past the valve seat. 

It is assumed that leakage past the valve seat is time dependent and the 
other two failure modes are cycle dependent. 

The probability of a single solenoid valve not failing closed in a 
single cycle of operation, RCSZ, is given by 

RCSZ = exp [ -LCLZ*10 -6 ] 

and the probability of no open failure or leak past the valve seat in a single 
cycle of operation, ROSZ^, is given by 

ROSZ fc = exp [ -LOPZ • 10 -6 - LVSZ • TIMEZ^ * 10 -6 ] 

LCLZ and LOPZ are the random or chance failure rates (failures per 10^ cycles) 
of a solenoid valve failing in a closed and open position, respectively. 

LVSZ is the failure rate (failures per 10 6 hours) of a solenoid valve due 
to excessive leakage past the valve seat. 

The solenoid valves may be arranged in a number of different configura- 
tions. A number of configurations are considered by the reliability model; 
namely the single valve ("S") , dual series ("DS") ,dual parallel ("DP”) , quad 
system ("Q") , and quad connected ("QC") systems. These configurations are 
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illustrated in Appendix III where reliability equations are derived in terms 
of the basic single valve reliability equations. Therefore, the reliability 
of the valve system is given by 
When VSC = "S" 

Then ' 

RSOLCZ, = ROSZ, + RCSZ - 1.0 
k k 

When VSC = "DS" 

Then 

RSOLCZ, = RCSZ 2 - [1.0 - ROSZ, ] 2 
k k 

When VSC = "DP" 

Then 

RSOLCZ k = ROSZ k 2 - [1.0 - RCSZ] 2 
When VSC = "Q" 

Then 

RSOLCZ k = [1.0 - (1.0 - ROSZ k 2 )] 2 - [1.0 - RCSZ 2 ] 2 
When VSC = "QC" 

Then 

RSOLCZ k = [1.0 - (1.0 - RCSZ) 2 ] 2 - [1.0 - ROSZ^] 2 

where VSC is the input variable which specifies the particular solenoid 
valve configuration to be used. It should be noted that RCSZ is not sub- 
scripted since only one cycle of operation is required. 

The reliability of the pressure transducer is assumed to depend 
primarily upon the time of actual use. The time of actual use, TIM^, is 
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the thrust time* and is computed by the performance portion of the model. 
Thus, the reliability of the pressure transducer is given by 
When HALFZ = 1 
Then 

RPDlZ k = exp [-LPZ 
RPD2Z k = 1.0 
When HALFZ * 2 
Then 

RPDlZ k - RPDlZ k _ 1 
RPD2Z k = exp [-LPZ 

where LPZ is the random or chance failure rate (failure/10 u hours) of the 
pressure transducer. RPDlZ k and RPD2Z k are the reliabilities of the pres- 
sure transducers in half-sys terns 1 and 2, respectively. 

It is assumed that the reliabiiity of the filter assembly, RFlZ k and 
RF2Z k for half-systems 1 and 2, respectively, are a function of both absolute 
time and use time. Absolute time must be considered because of the possibil- 
ity of failure of a welded connection. Therefore, the reliability of the 
filter assemblies is given by 


• • ip -6 ] 

3600 


where RPDIZ^ = 1.0 

• ^ • «r 6 ] 

3600 


*Since the performance model considers each half-system seperately, 

TIMj^ is the thrust time associated with the half-sys tern being evaluated. 

TIMj^ is initialized to zero after completion of the performance evaluation 


of the first half-system. 
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When 


HALFZ = 1 


Then 

RF1TZ = exp [-LFZ • TP V * 10 -6 ] 

3600 

RFlZ fe = RFlTZ k • exp [-LWZ • TIMEZ k • 10~ 6 ] 

{ 

RF2Z fc = exp [-LWZ • TIMEZ k • 10 -6 ] 

When HALFZ = 2 
Then 

RFlTZ k = RFlTZ k _ x where RFITZ^ =1.0 

RFlZ fc = RFlTZ k ' exp [-LWZ * TIMEZ fc • 10 -6 ] 

RF2Z = exp [-LFZ • TIM k • 10~ 6 - LWZ • TIMEZ. • 10~ 6 ] 

* 3600 k 

where LFZ is the chance failure rate (failures/10^ hours) of the filter media 

The performance portion of the model computes the number of impulses 

or cycles (J k ) required to achieve the velocity increment of the k th orbit 

correction. The total number of cycles of operation of half-systems 1 and 2 

(CYCLlZ k and CYCL2Z k , respectively) are thus given by 

When HALFZ =1 

When k = 1 

Then CYCL1Z. = J. 

k k 

CYCL2Z. = 0 
k 

When k > 1 

Then CYCL1Z, = CYCL1Z, _ + J. 

k k-1 k 

CYCL2Z k =0 
When HALFZ =2 
When k = 1 

Then CYCLlZ k =0 


CYCL2Z, = J, 
k k 
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When k > 1 


Then CYCL1Z. = CYCL1Z, , - 
k k-1 

CYCL2Z. = CYCL2Z. _• + J, 
k k-1 k 

The reliability equations for the thruster solenoid valve systems are 
similar to the half-system connecting valve. The single exception is that the 
thruster solenoid valves will experience many cycles of operation. The proba- 
bility of no closed failure of a single valve is 
RClZ fc = exp [-LCLZ • CYCLlZ k • 10 -6 ] 

RC2Z fc = exp [-LCLZ • CYCL2Z k • 10 -6 ] 

and the probability of no open failure or leak past the valve seat (of a 
single valve) is r 

R01Z k = exp [-LOPZ • CYCLlZ k • 10~ 6 - LVSZ * TIMEZ k • 10 — 6 ] 

R02Z fc = exp [-LOPZ • CYCL2Z k • 10~ 6 - LVSZ * TIMEZ k • 10 -6 ] 

where RC1Z AND R01Z refer to the first half-system and RC2Z and R02Z refer 

to the second half-system. As before, several different valve system types 
are allowed. These are specified by VST. . The valve system reliability is 
given by RST1Z for the first half-system and RST2Z for the second half- 

i 

system. Therefore (refer to Appendix III) , 

When VST = "S" 

Then 

RST1Z, = R01Z, + RC1Z, - 1.0 
k k k 

RST2Z. = R02Z + RC2Z - 1.0 
k k k 

When VST = "DS" 

Then 

RST1Z = RC1Z 2 - [1.0 - R01Z, ] 2 
k k k 

RST2Z k = RC2Z k 2 - [1.0 - R02Z k ] 2 
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When 


VST = "DP" 

Then 

RSTlZ k = R01Z k 2 - [1.0 - RClZ k ] 2 

RST2Z, - R02Z, 2 - [1.0 - RC2Z, ] 2 
k k k 

When VST = "Q" ' 

Then 

RSTlZ k = [1.0 - (1.0 - R01Z k ) 2 ] 2 - [1.0 - RClZ k 2 ] 2 

RST2Z k = [1.0 - (1.0 - R02Z k ) 2 ] 2 - [1.0 - RC2Z k 2 ] 2 

When VST = "QC" 

Then 

RSTlZ k = [1.0 - (1.0 - RClZ k ) 2 ] 2 - [1.0 - R01Z k 2 ] 2 

RST2Z fc = [1.0 - (1.0 - RC2Z k ) 2 ] 2 - [1.0 - R02Z fc 2 ] 2 

In a similar manner , the reliability of the pressurant and propellant 
tank solenoid valve system can be established. 

When VSP - "S" 

Then 

RSP1Z = R01Z. + RC1Z. -1.0 
k k k 

RSP2Z. = R02Z, + RC2Z. - 1.0 
k k k 

When VSP - "DS" 

Then 

RSP1Z, = RC1Z. 2 - [1.0 - R01Z, ] 2 
k k k 

RSP2Z fc = RC2Z k 2 - [1.0 - R02Z k ] 2 
When VSP = "DP" 

Then 

RSP1Z. = R01Z. 2 - [1.0 - RC1Z, ] 2 
k k k 

RSP2Z k = R02Z k 2 - [1.0 - RC2Z k J 2 
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When VSP = "Q" 

Then 

RSPlZ k = [1.0 - (1.0 - R01Z k ) 2 ] 2 - [1.0 - RClZ k 2 ] 2 

RSP2Z k = [1.0 - (1.0 - R02Z k ) 2 ] 2 - [1.0 - RC2Z k 2 ] 2 

. / 

When VSP = "QC" 

Then 

RSPlZ k = [1.0 - (1.0 - RClZ k ) 2 ] 2 - [1.0- R01Z k 2 ] 2 

RSP2Z k = [1.0 - (1.0 - RC2Z k ) 2 ] 2 - [1.0 - R02Z fc 2 ] 2 

In the above equations, RSPlZ k refers to the first half-system and RSP2Z k 
refers to the second half-system and VSP is the input variable which specifies 
the type of valve configuration to be considered. 

The thrusters contain a number of catalyst screens with failure 
rate L'SCZ. The reliability of the catalyst screens is a function of usage 
time. The thruster also contains a number of brazed joints (failure rate 
given by LBRZ) and weld joints. The reliability of the two thrusters is 
given by RTHlZ k and RTH2Z k . Therefore, 

When HALFZ =1 

Then , 

RTHZ = exp[-NSCZ • L'SCZ • TIM k • Iff 6 ] 

3600 

RTHlZ k = RTHZ k • exp [ (-NBZ • LBRZ * TIMEZ k -2 • LWZ- TIMEZ k ) • 10~ b ] 

RTH2Z k = exp [ (-NBZ • LBRZ • TIMEZ k - 2 • LWZ* TIMEZ k ) • 10 -6 ] 

When HALFZ =2 
Then 

RTHZ k = RTHZ k _ 1 where RTHZ^ =1.0 

RTHlZ k = RTHZ k * exp [ (-NBZ • LBRZ- TIMEZ k -2 -LWZ- TIMEZ k > • l<f 6 ] 

RTH2Z = exp [ (-NBZ LBRZ • TIMEZ -2 * LWZ- TIMEZ. - NSCZ 

• LSCZ • ^ • 10-6, k 
3600 
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The last component to be considered is the pressurant and propel- 
lant tank bladder. The function of the badder is to separate the pres- 
surant and the propellant. The bladder can fail due to the random or 
chance failure rate, LBZ, of holes developing in the bladder. The bladder 
is also subject to wearout phenomena since the bladder is forced to con- 
tract and expand for several cycles of operation. Prior to the need for 
repressurization of the first half-system, that is 
When M0DE1 <2 
Then 

RBlZ k = exp [-NPZ * LBZ • TIMEZ k • 10" 6 ] 

RB2Z. = RB1Z 
k k 


When repressurization of the first half-system is required, the reliability 
equations become 

When M0DE1 = 2 and M0DE2 <2 


Then 


RB1Z, 


= | exp [-NPZ * LBZ • TIMEZ * 10 6 ]} > 

. ' /*0O 

* < 2 J [exp(-(AZ - MBZ) 2 / (2.0 • STDZ 2 ))] 
/{J [exp(-(AZ - MBZ) 2 / (2.0 * STDZ 2 ))] 


dAZ 

dAz! 


RB2Z k = exp [-NPZ ' LBZ • TIMEZ k • 10" 6 ] 


where the ratio of the integrals indicates the probability of a second 
successful contract/expand cycle given that it survived the first. The 
first cycle is due to the initial filling of the tanks. When repressuriza- 
tion is required of the second half-system, the reliability equations are 
given by 



When M0DE1 = 2 and MODE 2 = 2 
Then 


RB1Z, 


= {exp [-NPZ * LBZ • TIMEZ • 10 _6 ]> 

2 J[ exp (-(AZ- MBZ) 2 / (2.0 • STDZ 2 ))] 
j^yiexp (-(AZ- MBZ) 2 / (2.0 * STDZ 2 ))] 


RB2Z. = RB1Z. 
k k 


dAZ 

dAZ 


AZ^ 


In the above equations it is assumed that wearout phenomena are 
adequately characterized by the normal distribution where MBZ is the ex- 


pected number of cycles and STDZ is the standard deviation of the number 
of cycles to failure. 

Before getting into the general system reliability equations, two 
configuration decisions must be made. A decision (SPFVZ) must be made as 
to whether the configuration will contain a single propellant fill valve 
for each half-system or a single valve per propellant tank. A decision 
(CONVZ) must also be made as to the type of valve system to be used for 
the half-system connecting valve. Therefore, 

^ ien SPFVZ = "Y" (i.e., a single propellant fill 

valve is to be used for each half 
Then system) 


RFVlZ k =1.0 

RFV2Z. = RFVZ. 
k k 

Whe* 1 SPFVZ = "N" (i.e., a single propellant fill 

valve is to be used for each pro- 
Then pellant tank) 

RFVlZ fc = (RFVZ^^ 2 

RFV2Z. =1.0 
k 
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Wh en CONVZ = "SQ" or not "SOL" (i.e., specifications of the 

type of valve system connecting 
Th en the two half-systems) 

BZ, = RSV2Z, 
k k 

When CONVZ = "SOL" 

Then 

BZ. = RSOLCZ, 
k k 

The overall system reliability can now be determined in terms of the 
reliability of the various components and the specific arrangement of the 
components. In general three basic situations are considered: 

(a) repressurization is not required for either half-system 

(b) repressurization is required for the first half-system, 

(c) repressurization is required for both half-systems. 

The probability of having a working thruster system available for use when 
using the first half-system is 

RElZ fc = 1.0- jl.0 - RFV2Z k • RPDlZ k * RFlZ fc • RSTlZ k • RTHlZ fc • RLTZ k | 

1.0 - BZ • RFV2Z, • RPD2Z. • RF2Z. -RST2Z, • RTH2Z, • RLTZ, 
k k k k k k k 

and represents the probability of thruster system #1, or thruster system #2 
and half-system interconnecting valve functioning properly. Similarly, the 
probability of having a working thruster system available for use when using 
the second half-system is * ' 3 

RE2Z fc = 1.0 - | 1.0 - BZ fc ’ RFV2Z k ‘ RPDlZ fc • RFlZ fc • RSTlZ k • RTHIZ^RL^ 

• | 1.0 - RFV2Z k • RPD2Z k • RF2Z k • RST2Z 
The probability of successful operation of a propellant tank and its 
associated fill valve, solenoid valve, and lines is given by R1Z for the 

(V 

first half-system and R2Z fc for the second half-system. Therefore, 


, -RTH2Z. 
k k 


RLTZ, 



28 



R1Z, = RB1Z • RFV1Z, • RSP1Z, • RLPPZ, 
k k k k k 

R2Z. = RB2Z, • RFV2Z. • RSP2Z, • RLPPZ, 
k k . k k k 

QlZ k = 1.0 - RlZ fc 
Q2Z fc = 1.0 - R2Z k 

The reliability of the high pressure tank, fill valve, and squib 


valve system is given by DELlZ k and DEL2Z k and depends upon the need for 
repressurization. Therefore, 

When M0DE1 < 2 


Then 


DELlZ k =1.0 
DEL2Z k =1.0 

When MODEl = 2 and MODE 2 < 2 

Then 


When 


DEL1Z. = RFVZ, 
k k 

DEL2Z fc =1.0 
MODEl = 2 and 


• RTZ • RSV1Z . 
k 

MODE 2 = 2 


Then 

DELlZ k = RFVZ k • RTZ k • RSV1Z ... 

DEL2Z. = DEL1Z, 
k k 

Since there are multiple propellant and pressurant tanks and redundant 
thruster systems, the reliability of. the system is a function of the propel- 
lant requirement as a function of time. The model computes the propellant 
required for each maneuver and thence establishes the reliability in terms 


of the probability of the propellant being available. The possible ways of 
achieving six different outcomes, or propellant levels, are considered. 
Specifically, the probability of having 1/6, 1/3, 1/2, 2/3, 5/6, and full 
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propellant capability at any instant of time is determined. The way of 
achieving these outcomes is illustrated in Figures 2-7 where the following 
nomenclature applies: 

A & B - refers to half-system one and two. 

1A, & IB - refers to high pressure tank system, fill valve, 
and squib valve system. 

. 2A, & 2B - refers to fill valve, lines between squib and 

propellant tanks, and propellant tanks. 

. 3A, 4A, 5A, 3B, 4B, 5B - refers to the propellant tank 

bladder, fill valve (A), solenoid valve system, 
and lines between propellant tanks and fill valve (B) 
There are six tank systems; three per half. 

. 6 - refers to the half-system connecting valve system. 

. 7A, 7B - refers to fill* valve (B) , pressure transducer, fil- 

ter, orifice, solenoid valve system, and thruster. 
Using figures 2-7 as a guide, the probability of having different levels 
of propellant available when needed and functioning properly can be estab- 
lished as 

GZ, = RLHPZ, • RFVZ, • RPTZ, * 

k k k k 

PFZ fc ± = ^3.0 * GZ fc | ' |klZ k • QlZ k 2 • RElZ k • DELlZ k • [1.0 - GZ fc • DEL2Z k 
+ GZ fc • DEL2Z k • Q2Z k 3 ] + R2Z k • Q2Z^ • RE2Z k • [1.0 - GZ fc • DELlZ k 

+ GZ k • DELlZ k • QlZ k 3 ] } _ * 

PFZ k,2 = I 3,0 ' GZ k } ' { R1Z k 2 ’ Q1Z k * RE1Z k * DELlz k * El-0 - GZ fc • DEL2Z k 
+ GZ k * DEL2Z k • Q2Z k 3 ] + R2Z k 2 • Q2Z k • RE2Z fc • DEL2Z k • [1.0 - GZ fc 

• DELlZ k + GZ k • DELlZ k QlZ k 3 ] + 3.0 • RlZ fc • R2Z fc * QIZ^ 2 

• Q2Z, 2 • RE1Z, • RE2Z • DEL1Z, • DEL2Z, • GZ, \ 

k k^k k k k) 
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FIGURE 5 
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PFZ k,r“ { GZ k}’ { R1Z k 3 ' RE1Z k * DEL1Z k * [1<0 " GZ k ’ DEL2z k + GZ k • DEL 2 Z k 

. Q 2z k 3 3 + R2Z k 3 * R E2Z fc • DEL2Z k • [1.0 - GZ k • DELlZ k + GZ k • DELlZ fc 
• QlZ k 3 ] +9.0 • RlZ k • R2Z k • Q2Z k • RElZ fc * RE2Z fc * DELlZ k * QlZ fc 
° DEL2Z k • GZ k • [ Rlz k * Q 2z k + Qiz k • R2Z k ] | 

p EZ k ^ => ^3.0 RElZ k •, RE2Z fc • GZ k 2 • DELlZ fc • DEL2Z k • RlZ fc • R2Z fc j 

•| Rlz k 2 * Q 2z k 2 + Glz k 2- R2Z k 2 + 3.0 ’ RlZ k • QlZ fc • Q2Z fc • R2Z fc J 

PFZ k>5 = ^3.0 RElZ k • RE2Z fc • GZ^. 2 • DELlZ k • DEL2Z fc • RlZ fc 2 • R2Z fc 2 J- 

•{ R1 \ • ^ z k + « lz k ' R:!Z u} 

PFZ fc ^ 6 = RElZ k • RE2Z k • GZ fc 2 • DELlZ k • DEL2Z k • RlZ fc 3 • R2Z fc 3 

where pp Z k is the probability of the system functioning properly and having 
exactly i propellant tanks available. Therefore, the probability of having 
more than the desired propellant available when needed and the system func- 
tioning properly is given by p MZ k . 

When HALFZ =1 


(i.e., at least 1 propellant tank 
available) 


Then PMZ fc = PMl k + PM2 k 
When HALFZ = 2 


Then PMZ fc = PM2 fc 

When PMZ > 5.0 • PMO * 2.0 . 

6.0 


Then PZ, = 

k 

IZ=1 



— H ■ FM° 
6.0 

Then PZ, = 
k 

• PMO 

3.0 

Then PZ, = 

k 


2.0 > PMZ > 1A. 

3.0 


PMO * 2.0 (i.e., at least 2 pro- 

pellant tanks available) 


t 


IZ=2 


PFZ 


k,IZ 


2.0 > PMZ > 1.0 '■» PMO • 2.0 

2.0 

6 


(i_.e^ r at least 2 pro- 
pellant tanks available) 


PFZ 


k,IZ 


IZ=3 
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-^4 • PMO • 2.0 > PMZ > _lj_0 • PMO '• 2.0 

2.0 6 k 3.0 

ass. PZ k " 2 PFZ k,iz 

IZ=4 


(i.e., at least 4 pro- 
pellant tanks available) 


When 1.0 • PMO • 2.0 > PMZ > 1^0 

3.0 6.0 

6 

Then PZ k = PFZ k,IZ 

IZ=5 


• PMO * 2.0 (i.e., at least 5 pro- 

pellant tanks available) 


When 1.0 • PMO • 2.0 > PMZ, > 0 (i.e., all 6 pro- 

6.0 pellant tanks available) 

Then PZ t = PPZ t>6 

In other words, PMZ^ is the reliability of the orbit correction system as 
a function of the orbit correction maneuver. The overall system reliability 
takes into account the specific operations of the auxiliary propulsion system; 
i.e., thrust duration, number of impulses to achieve desired velocity incre- 
ments, magnitude and timing of velocity increments, etc. The reliability 
model considers each of the k maneuvers as a different outcome and establishes 
the probability of achieving each outcome taking into account the various ways 
in which the outcome can be achieved. 

Typical Results 

In order to demonstrate the type of results available from this model, 
a typical mission profile was assumed along with propulsion system character- 
istics. The set of input data is printed out and is attached as part of Ap- 
pendix II. The mission evaluated consisted of 39 orbit corrections spaced 
over a period of approximately one year . Several of the orbit corrections are 
spaced at intervals of one-half an orbit period while others are seperated by 
tens of days. The computed results are also shown in Appendix II. Only 37 
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maneuvers are possible due to the performance of the propulsion system. 
The results of the reliability computations are summarized in Figure 8 
which indicates the probability of failing to perform the desired man- 
euver in terms of the maneuver number. This is one minus the probability 
of success or the "unreliability" of the system. The specific numbers 
are unimportant. What is important is the rapid manner in which the 
probability of a failure increases as the number of maneuvers increases. 
The reason for this is that early in the mission only one or more pro- 
pellant tanks need function successfully. In the latter stages of the 
mission, maneuvers can only be performed if all six propellant tanks 
have and continue to function properlyi The large propability incre- 
ments (k=21, 25, 32) are due to the requirement of additional propellant 
tanks being available. The small changes in the probability of failure 
between maneuvers no. 25 & 26, 29 & 30, etc., are due to the fact that 
the maneuvers occur within a time interval of 1/2 an orbit period (i.e., 
a perigee followed by an apogee correction) . 

GENERAL DISCUSSION 

A general methodology for evaluating the reliability of a propulsion 
system as a function of time or mission objective has been described. 

The method requires the interaction of a propulsion system performance 
model with a reliability model. The performance model provides basic 
timing and cycle data required for the reliability computations and the 
reliability model provides the basic structure for the manner in which 
the various mission outcomes or objectives can be achieved. The relia- 
bility model thence performs the reliability computations, in terms of 
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specific equipment configurations, providing probability of success data 
for achieving each of the desired mission objectives. 

When a mission consists of multiple time distributed mission objectives 
(for example, the previously discussed orbit correction system), the nor- 
mal procedure for evaluating system reliability at a single point in time 
loses significance. The reliability of the system after t hours is still 
an important indicator of overall system performance. However, it is also 
important to establish the probability of successfully achieving each de- 
sired mission objective within the t hours. 

The availability of data indicating the probability of successfully 
achieving each of the desired mission objectives can provide valuable in- 
sights for the comparison of system alternatives. In particular, if the 
relative importance of the various mission objectives is established, a 
figure of merit can be computed and used to assist with the ranking and 
comparison of alternatives. The figure of merit is of the form 

F = • I.M?. 

k 

where IMP^ is the relative importance of achieving the k fc ^ mission objec- 
tive. It is normally desirable to chose that system configuration which 
maximizes F . When cost constraints become important, the above can be 
utilized to establish a measure of cost effectiveness (CE) such that 

CE - F 
C 

where C is the total cost associated with the system configuration under 
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study. The objective is to chose that system configuration which maximizes 
the value of CE. When costs are spread over a significant period of time, 
present value concepts should be used so as to take into account the timing 
of expenditure of funds. _ - 
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APPENDIX A: Program Listing 


A-l 


HEAL LVZ,LCZ,NTZ,LWZ,NWHPZ,LCLZ, LOPZ, LVSZ,LPZ ,LFZ, 

1NS CZ.LSCZ, N BZ , LBR Z , NP Z , LEZ , NWTZ , NW PZ _ 

REAL*8 MBZ , STDZ , RT2S, ERFC2, ERFC1 , ERFC21 , DERFC 
INTEGER CYCL1Z, CYCL2Z 

IHTJE.GER._HA L FZ _____ 

DATAKS/ * S */, KDS/ ' DS ' / , KDP/ ' DP */ , KQ/ * Q* / ,KN/* N */, K SQ/* SQ */ 

DIMENSION ORF (10) , DELE (10) , DEL PE (10) f DEL V (100) , CD (10) ,A (10) , 

„1STI(.100) ,STIME (100) ,SDELV.(100) ,STDEL (100) ,STIM(100) 1 _ 

DIMENSION SFVZ (100) / RTZ(100) , RPTZ (100) , RLBPZ (100) , 

1 TIMEZ (10 0) , RF1TZ (100) ,RPD1Z (100) ,RPE2Z (100) ,PFZ( 100,6), 

2RF1Z (100) , RF2Z (100) , CYCLlZ(IOO) ,CYCL2Z (100) ,PZ(100)., 

3RE1Z (100) , RS2Z ( 100 ) , RSP1 Z ( 1 00) ,RST1Z (100) ,RST2Z (100) , 

4RSP2Z (100) , RTH1 Z (100) , RTH2Z(1Q0) , RB 1 Z ( 1 00 ) , RB 2Z ( 1 00) , BZ ( 1 00) , 

5RLPPZ (100) , RLTZ (100) ,RTHZ (100) .. _ __ 

NAME LIST/I NDATL/VE JIN, Z,CFI N , DEN ,AT , PFIL , STG , CD, A,DT ,DTI,P ,DPT, 

1 PMO , VIK, VTINT, PW,WSC,L,ULV,D£LV,PRES ,PCMIN,PCIN, DP, T MAX , TLNCH , 

2B GAS.G _ .I 

NAME LI ST/I NDATG/TIMEZ ,LVZ,LCZ , NTZ , LWZ , NS 1 Z , PFS Z , 

1NWH?Z,C0NVZ,NS2Z ,LCLZ ,LOPZ , LVSZ ,VSC ,LPZ,LFZ ,NSCZ, 

2LSCZ , NBZ , L3RZ ,NPZ , LBZ , 8 BZ , STDZ , N WTZ , N.WPZ, SPF VZ ,YSC, VST,.VSP : 

WHITE (6 ,70) 

70 FORMAT (• 1INPOT DATA 1 ) 

C FOB STEADY STATE MODE OSE PW= 1 . OE+ 1 0 .. SEC. ' - 

READ (5, INDATL) 

WRITE(6,INDATL) 

_READ (5, INCAT G) . 

WRITE (6, INDATG) 

C AT INPOT DIVIDE THE FOLLOWING VARIABLES BY 2 

- P MO= PMO/2. ’ 1 1 ; 

VTK= VTK/2 . 

ULV= OL V/2. 

V T INT=VTINT/2. . .. - ■ - J 

C COMPLEMENT OF THE ERROR FUNCTION 
HT2 S=1 .4 1 4 13562 3731D0*STDZ 

E RFC2 = DERFC( (2 . DO-MBZ) /RT2S) 

ERPC1=DERFC( (1 . DO-MBZ) /RT2S) 

ERFC21=ERFC2/ERFC1 » 

C ALL . STINT (0 . , 0 . , 0 . , 0 . , - 1 , N GR IP E, 0 , 0) ... 1 1 

IF (NGRIPE. NE.Q) GO TO 100 
WRITE (6, 85) 

86 F ORMAT (V BNP DATA' ,///) ; — 

WRITE (6, 83) 

83 FORM AT ( * LIST OF VARIABLES WITH CORRESPONDING UNITS*,/) 

WRITE (6,82) 

82 FORM AT (' AT=THRO AT AREA (SQ MET) PGO=INIT TANK PRES(N/SQ MET) STG-S 
IT AGE NUMBER CD=ORIFICE COEF A=ORIFICE OR VALVE FLOW AREA (SQ MET) ' 
2 /* WSC=SP ACECR AFT MASS... (KG) U LV =GAS. BOT...VOLUME _(C U M ET) DEL V= CELT A .... 
3V (MET/SEC) ') 

WRITE (6,84) 

_84 FORMAT ( * Z=CCMPRES FACT CF=THRUST COEF DEN = PROP DENS (K/CUMET) 

1 DT=LOOP TIME INCH (SEC ) DT I=LOO P2 TIME INC (SEC) P=PRES BOT EXPAN 
2EXP */ * DPT=I NTERN AL TANK PRES DROP (N/SQ MET) PMO=PRCP MASS (KG) 

3R=PROP TANK RADIUS (MET) ... PW= PULSE .WIDTH (SEC) ') 

WRITE (6,81) 

81 FORMAT (' VTK=TANK VOL (CU MET) VT INT=INTERNAL TANK DISPLACEMENT (CU 
1 MET) VG2 +_ VG 3 =.. 2N D . ST 6 G AS VOL. (CUMET) A(2)=ORIFICE AREA (CUMST 



2) •) 

PMI=PMO 

C*** THIS ROUTINE ADJUSTS THE TANK PRESSURE FOR LAUNCH TEMPERATURE 
TEMP =TMAX 

£_ALL_S TIN T_(T EM P ,0. f 0. P DEN 1 , NGRIPE ,5 f 5) ___ 

IF (NGRIPE. NE. 0) GO TO 100 
DENO=DEN 

VOL G = V.T K - V T.I NT-PMO/D E NO 

DEN=0. 

TEMP=TLNCH 

CAL L_S T I NT (T E M P ,_0 0., DEN, 1, NGRIPE, .5 ,3) 

IF ( N GRIPE . NE . 0) GO TO 100 
DELV CL= (PMO/DENO) * (1 . -DENO/DEN) 

V G 0 = V OLG+DELVOL ' 

PGO=PFIL*VOLG*TLNCH/ (TMAX+VGC) 

C2=Z/DEN 

C._ TiHS_ROUT IN E_S IZ ES_T_HE_EAO.N__.OR IF_I.CE ' 

V EJ=VEJIN 
. CF=CFIN 

AP=_._7_8 53975* DP* DP_ ______ . 

PG = PGO 

DPT=DPT*6897.4862 

D O 131 M= 1,5000 

SUMA=0. 

DO 10 1=1,10 

I_£_(_CD (I) . EQ. 0. ) GO TO l0 1 1 1 

IF (A (I) . EQ . 0 . ) GO TO 101 
ORF(I)=CD(I)*'CD(I)*A(I)’» ! A(I) 

JL0_SUHA=S UM A + ( 1-. -0 R r. ( I ) /. ( A P * A P ) ) /(O R F( I )_*2 ,_*EE N ) 

101 C3=SUMA 

KM=I-1 

CS.TAR =V EJ/C F 1 

C1=CSTAR/AT 

C0NST=(C1*C1) -4. *C 3* ( DPT-PG) 

FM=.(-C1 + SQRT (CONST) ) / (2. *C3) _____ 

IF (FE.LE.O.)GC TO 114 
PC=FM*Cl 

C ALL STINT (PC, 0. ,0. , VE J , 1 , NG RI PE , 2 , 2) 

IF (NGRIPE. NE. 0) GO TO 100 

CALL STINT (PC,0. ,0.,CF, 1 , NGRI PE , 3 , 3) 

_____I_F4.N.GR I PE. NE. 0) GO TO _1 00 

IF (PC.GE.PCIN) GO TO 133 
A (1 ) =A ( 1) +.5E-08 

IF (A (1) .LE. AP) GO_TO .131 

A (1 ) =AP 
GO TO 133 

1_3J_CO N T I N U E __ 

133 DOR=SQRT (A (1) *1.274324) 

DO 129 N = 1 , KM 

DELP (N) = (FM*FM) *(1 .-ORF(N) / (AP*AP) ) / (ORF (N) *2..*DEN) . 

129 DELPE(N) =DELP(N)/6897. 4862 
KI = 1 

HALFZ=1 , 

MODE1=0 

MODE2=0 

2 VEJ= VEJIN : 
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IERR=0 

CF=C F IN I — 

TIM=0. 

SUM= 0. 

ICTR=0 ; 

TDEL=0. 

TI=0 . i 

MODE=0 ; — ] 

RESTIM=0. | 

PM=PMO ! 

£**jL_L__IJS T H E .. N U M B E R__0 F_.M AN UE VERS : 1 

DO 40 K=KI , I ' 

TIMEP=0. J 

_ B.TIM E=R ESTIil ! : 1 — 

RESTIM=0. 

DELT AV=0 . 

DELN=0. : i ■: 

C THIS LOOP IS ONE POLSE , 

DO 30 J=1 ,10000 i 

__ BTIME=BTIME+TIMEP ._ ; j 

C THIS IS THE FIRST STAGE OF BLOWDOWN 

DO 20 1=1, 10000 | 

IF (MODE. EQ. 2) GO TO. 115. -| 

IF (HALFZ. EQ.2) GO TO 112 | 

MODE 1=0 

MODE 2=0. : ; 

113 IF (STG. EQ. 1 .) GO TO 15 i 

IF (PG-PRES) 11,11,15 } 

11 YGO=7G O+S 0 M*C2+OL V 

PGO^PFIL 

SOM=0. 

P HO= FM ' __ ; 

MODE=2 

C THIS IS BOTH STAGES OF BLOWDOWN 

15 CS TAR = V E J/CF_; 

Cl =CSTAR/AT ; 

C THE EXPONENT WAS EXPERIMENTALLY DETERMINED , | 

P_G= P GO * (VGO/ (VG0 + SUM*C2) ) **P 1 1 < 

PCT= (PMI-PM) /PMI I 

C THE TANK INTERNAL PRESSURE GRADIENT IS A FUNCTION OF REMAINING PROPELLANT ; 

C ALL STI NT (PCT, 0 . , 0 . , DPT , 1 , N GR IPE ,.1 ,1) I. L 1 : 

DPT=CPT*6897 .486 ' 

IF (NGRIPE.NE.O) GO TO 100 

C ON ST= (C1*C1)-4.*C3»(DPT-PG) 

IF (CONST. LE. 0.) WRITE (6,65) | 

65 FORMAT ('LOOK UP ERROR') ♦ 

FM= (-C1 +SQRT (CONST) ) /(2.*C3) • . J 

I F (FM. LE. 0. ) GO TO 21 j 

SUM=SUM+DT*FM 

■ _ PM=PMO- SUM - j 

IF (PM.LE.O.)GO TO 21 I 

TIME=DT*I 

T I M E P = T I M E 1 i 

PC=FM*C 1 | 

IF(PC.LE.PCMIN) GO TO 21 j 

_C THIS ..CURV E__V ARIES^FQR DIFFERENT .TH RUSTER S , , . ; J 



CALL STINT(PC,0. , 0 . , CF , 1 , N GRIP E , 3 , 3) 

IF.(NGRIPE. M B. 0) GO TO. 100 . . 

CALL STINT (PC, 0. ,0. , VEJ, 1,NGRIPE, 2,2) 

IF (NGSIPE.NE.O) GO TO 100 

C THIS C URV E _y A R I ES_.FO R DIFFER ENT. P_ U LS E_MQ.EES ' . _ 

IF(PW.EQ. 1.E+10) GO TO 8 
PN=J 

CALL STINT (PN, 0. ,0. , EFF, 1 , NGF.IPE, 4,4 ) . 

IF (NGRIPE.NE.O) GO TO 100 

VEJ=VEJ*EF? , 

8_8IT=FM* V E J * D T. ! 

TI=TI+BIT 
THR=FM* VEJ 

S_PI_= V E J / .9 . 8 _ . j 

IF (PW.NE. 1. E+10) GO TO 9 

- DELTAV=DELTAV+VEJ*ALOG( (WSC+PM+FM*DT) / (WSC+PM) ) 

IF.(DELTAV. GE. DELV (K) ) GO TO 31 ; 

9 IF (MODE. EQ. 2) GO TO 111 
C THIS FOR FIRST BLOWDOWN 

IF. (K.N E . 1.) GO TO 111 ; 

IF(J.NE. 1) GO TO 111 
IF (I.NE.1) GO TO 11 1 

£.O.R=FH*FH* (.1 . -ORE (1) / ( AP * AP ) ) / (OR F (1)*2.*EEN) 

PINLET=PG-DPT-POR 
WRITE (6,69) 

6.9._FQRM AT (*_1 INITIAL CONDITIONS * ,//) ... . 

WRITE (6,59) A ( 1) 

59 FORMAT (' ORIFICE AREA = •,E11.4,»'SQ MET*,/) 

W RIT.EL ( 6 , 6 2 ) DO R , 

62 FORM AT (* ORIFICE DIAMETER = *,E11.4,' SQ MET',/) 

WRITE (6, 64) PRES 

64. FORM AT (* .PRESSURE AT REPRESSURIZATION , E 1 2. 4,.'_N/SQ. MET*,/) 

IF(A (2) .NE. 0.) WRITE (6,60) DELPE(I) 

60 FORMAT (* ORIFICE PRES DROP = ' ,F5.0,* PSI',/) 

IF.(A (2) .NE. 0.) WRITE ( 6, 58) PINLET. 

58 FORMAT (• PINLET = *,E12.4,' N/5Q M * ,/) 

WRITE (6,61) j 

£-1_F_0 RM.AT ( !_M AN E U V ER * ,6X,.' PROP* , 12X, . 'TANK*., 12X, 'CHAM* , 14X,. . 

1 'THRUST*, 16X, 'VEJ* ,15X, *ISP*,11X, ' FLGW' 

2,/,* NUMBER', 7X, *LEFT',12X, *PRES*,12X, * PRES ', 68X, • RATE 

3.',/, 15X, • (KG) * ,10X, ' (N/SQ M) ' , 8 X ,*( N/S Q M )*_,14 X ,.* ( N ) .' ,15 X ,*( M/ S EC ) 

4' ,12X,' (SEC) ' ,3X, • (KG/SEC) '/) 

WRITE (6,71) K,PM,PG,PC, THR , V E J , SP I , F M 

7.1.-.. FORM AT (4X, 13 ,4X , 3 (E 12.4, 4X) ,E17. 4,4X , El 7. 5, 4X, F 1 0.2, 4 X, El 2 . 4 ,/) 

111 IF (I + DT.GE.PW) GO TO 1 9 
GO TO 20 

1J_2_MODE 1 =2_„ ; 

MODE 2=0 
GO TO 113 

1.1 5_I F (H AL FZ . EQ . 2) GO TO _1 1 6 ■ __ 

MODE 1 = 2 
MODE2=0 

GO. .TO.. 15 , ... : 

116 MODE 1=2 
MODE 2= 2 

GO TO . 15 
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20 CONTINUE 

19.. PBIT=SUM*VEJ ... . 

DELTAV=DELTAV+VEJ*ALOG ( ( WSC + PM+FW*DT) / (WSC+Pti) ) 

IF (DELTAV. GE. D3LV (K) ) GO TO 31 

10... CONTINUE .. i 

RES1IM= BTIKE+TIM EP 7 

31 NJ=J 

IF (HALFZ . EQ. 2) GO TO 1 000 . ... : : 

C NUMBER OF CYCLES OF THE DELTA 7 THRUSTERS ; HALFZ=1 
IF (K. EQ.1 ) GO TO 1010 

CYCL1Z (K) =CYCL 1 Z (K-1 ) + NJ 

CYCL2Z (K) =0 
GO TO 1064 

1 01 0 CY CL1Z (K1 =N J ■ 

CYCL2Z (K) =0 
GO TO 1064 

C NUMBER .OF CYCLES OF T H E..D E L T A. JLTH.RU S TERS; .H ALF.Z= 2 

1000 IF (K . EQ. 1) GO TO 1016 ; 

CYCL1Z (K) =CYCL1Z (K-1) 

C.YCL2Z (K) =CYCL2Z (K-Jl +N J 

GO TO 1064 . 

1016 C YCL 1Z (K) = 0 

C YCL2Z (K) = NJ . : 

1064 TDEL=TDEL+DEL'T A V 

POR=FM*FM* (1.-ORF (1)/ (AP + AP) )/ (ORF(1) *2.*DEN) 

PIN LET=PG-DPT-POR 1 

IF (ICTR. EQ. 0) GO TO 16 
IE (ICTR. GE. 54) GO TO 16 

6 W R I T E (.6,71) K, PM, P G , P C , I H R ,VE J , S PI ,.F. M 

ICTR=ICTH+2 

C SAVE THE FOLLOWING VARIABLES FOR LATER PRINTING 
S.TI.(K) =T.I 

IF (NJ. GT. 1) STIME (K) =BTIME 
IF (NJ.EQ.I)STIMS (K)=TIME 

SDELV (K) =DELTAV ; ! i ; 

STDEL (K) =TDEL 
IF (NJ.GT. 1) TIME=BTIME 

_TJM_=.TIM+TIME , .... * 

STIM (K) =TIM 
NK = K 

: TIME Z (K) =24. *TI MEZ (K) .J 

C RELIABILITY OF FILL VALVES 

1 001 RFVZ (K) = 1 (1.-EXP (-LVZ*TIMEZ (K) *1. OE-6) ) * (1.-EXP (-LCZ*TIMEZ (K) 

1*1.01-06)) ... ' 

C RELIABILITY OF HIGH PRESSURE TANK SYSTEM 
RTZ (K)=EXP(-NTZ*LWZ*TIMEZ (K) *1 .0E-06) 

_C_RE.LI ABILITY OF SQUIE VALVE SYSTEM (SV1) ; ' 

RSV1 Z=1 (1 .-PFSZ) »^NS1Z 

C PROPELLANT TANK WELD RING RELIABILITY 

RPTZ (K)=EXP (-LWZ+TIMEZ (K) *1 .0E-06) 

C RELIABILITY OF WELDED CONNECTIONS BETWEEN SQUID VALVE (S V 1 ) 

C AND PRESSURANT AND PROPELLANT TANKS 

RLH.PZ (K). = EXP (-N WHP Z^L WZ^TTM EZ (K) *1 . 0S-06) 

C RELIABILITY OF SQUIB VALVE SYSTEM (SV 2) 

IF (CCNVZ. NE. KSQ) GO TO 1002 

RSV2 Z=J..-(1.-P FSZ)**NS2Z 
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GO TO 1003 

_C_HJE LIABILITY. OF SOLENOID VALVE SYSTEM (S.CLC) 

1002 RCSZ=EXP (-LCLZ*1 .0E-06) 

ROSZ=EXP (-LOPZ* 1 . 0E-06 -LVSZ*TIMEZ (K) *1 . 0E-06) 

IF (7 SC. EQ . KS ) GO TO 1 004 _ 

IF(VSC. EQ.KDS) GO TO 1005 
IF (VSC. EQ . KDP) GO TO 1006 

I F (VSC. EO.KO) GO TO 1007 ' 1 

G QUAD CONNECTED VALVE SYSTEM 

RSOLCZ= (1 .-(1 .-3CSZ) * (1.-RCSZ) ) **2- (1 ,-ROSZ*ROSZ) **2 

GO ... TO. 1 0 0 a 

C SINGLE VALVE SYSTEM 
1004 RSOLCZ= ROSZ+RCSZ- 1 . 

G O TO 1 008 

C DUEL SERIES VALVE SYSTEM 

.1005 RSOLCZ=RCSZ*RCSZ- (1.-ROSZ) * (1-ROSZ) 

G..Q_T.O 1 008 

C DUEL PARALLEL VALVE SYSTEM 

1006 RSOLCZ=ROSZ*ROSZ-( 1 ,-ECSZ) * (1 .-RCSZ) 

G O T O._10 0 9 

C QUAD VALVE SYSTEM 

1007 RSOLCZ= (1.- ( 1.-ROSZ) * (1.-R0SZ) ) **2- ( 1 . -RCS Z*RCSZ) **2 

_C_B EL LABILITY 0 F .PRESSURE. TRANSDUCE S 

1008 IF (HALFZ. EQ. 2) GO TO 1009 
C HALFZ EQUALS ONE 

_ RPD 1 Z ( K) = EX P(-LPZ* (TIM/3 600.) * 0E- 06 ) 

RPD2Z(K) = 1. 

C RELIABILITY OF FILTER ASSEMBLY; HALFZ=1 

. R F1TZ (K) = EXP (LFZ* (TIM/3600. ). * 1. 0E-O6) 

RF1Z (K) =RF1TZ (K) *EXP ( -LWZ* II MEZ (K)*1.0E-06) 

RF2Z (K) =EXP(-LWZ*TIMEZ (K) *1. 0E-06) 

GO _ TC .1 0 11 ; ^ __r 1 

C HALFZ EQUALS TWO 

1009 IF (K . EQ . 1 ) GO TO 1012 

i HPD.1.Z (K) =RPD1Z (K-L) 

GO TO 1013 
1012 RPD1Z (K) =1 . 

_.UQ.13_R PD 2 Z ( K) = EXP (- L P Z* ( T IM/ 3 6 0 0 . ) * 1 . 0 E- 0 6 ) * 

C RELIABILITY OF FILTER ASSEMBLY; HALFZ=2 
I F (K. EQ. 1) GO TO 1014 

R F1IZ (K) =RF1 TZ (K-1 ) 

GO TC 1015 

1014 RF1 TZ (K) =1 . „ 

__1. 01.5, RFlZ (K) =RFlTZ (K) *EXP ( -LWZ*TIMEZ (K)*1 .0E-06) 

RF2Z (K) = EXP (-LFZ* (TIM/3600. ) *1. 0E-06 -LWZ *TI MEZ ( K) *1. 0E-06) 
C RELIABILITY OF THRUSTER SOLENOID VALVE SYSTEM 

_C PROBABILITY... OF NO CLOSED FAILURE OF _.A_SINGLE._VALV.E 

1011 R C 1 Z = EX P (-LCLZ *CYCL1 Z ( K) *1 . 0E-06) 

RC2Z=EXP (-LCLZ*CYCL2Z (K) *1. OE-06) 

C PROBABILITY CF NO OPEN FAILURE OR LEAK PAST VALVE SEAT 

R01Z=EXP(-L0PZ*CYCL1Z (K) *1. 0E-06-LVSZ*TIM EZ (K) *1 . 0E-0 6) 
R02Z=EXP (-LOPZ *C YCL2Z (K)*1. 0E-06 -LVSZ*TIMEZ(K)*1. 0E-06) 

_XF ( VS T . E Q . K S ) GO TO 1017 . L _ 

IF (VST. EQ. KDS) GO TO 1018 
IF (VST. EQ.KDP) GO TO 1019 

IF .( V S I ... EQ .. K Q )GO T0...1 0 20 ; 
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C QUAD CONNECTED VALVE SYSTEM 

g ST 1 Z_(K ) = ( 1 . (1. - RC 1 Z ) * { 1 . -R C 1 Z) ) ** 2- ( 1 • - S0 1 Z *F0 1 Z) **2 _ 

HST2Z (K) = (1.- (1 . -RC2Z) * (1.-RC2Z) ) **2 - ( 1 . -RO 2Z*RO 2Z) **2 
GO TC 1021 

__C SINGLE VALVE SYSTEM . ; - 

1017 RST1Z (K)=R01Z + RC1Z-1. 

RST2Z (K) =R02Z+RC2Z-1 . 

G O TO 1021 _ . _ 

C DUEL SERIES VALVE SYSTEM 

1018 RST1Z(K)=RC1Z*RC1Z-(1.-R01Z)*(1.-H01Z) 

RST2.Z(K) =RC2Z*RC2Z- (1.-R0 2Z) * (1.-R0 2Z) 

GO TO 1021 

C DUEL PARALLEL VALVE SYSTEM 

1019 R ST 1 Z (K) =R0 1 Z*RO 1 Z- (1.-RC1Z) * (1.-RC1Z) __ ■ 

RST2Z (K) =R02Z«R02Z-(1 .-RC2Z) * (1 .-RC2Z) 

GO TO 1021 

__C__aUAD.. VALVE SYSTEM ; 

1020 RST1Z (K) = ( 1. -(1.-R01Z)* (1.-R01Z) ) ** 2 - <1.-RC1Z*RC1 Z) **2 
RST2Z (K)= (1.- (1 .-R02Z) * (1.-R02Z) ) ** 2 - ( 1 . - RC 2Z*RC2Z) **2 

C TH RUSTER RELIABILITY _ _ 

1021 IF (HALFZ. EQ. 2) GO TO 1022 

RTHZ (K) =EXP (-NSCZ*LSCZ* (TIM/3 600.)*1.0E-06) 

RIH2Z_(K)=EX? ( (-NBZ*LBRZ*TIMEZ (K)-2. *LWZ*TIMEZ (K) ) *1. 0E-06) 

RTH1Z (K) =RTHZ (K) *RTH2Z (K) 

GO TO 1023 

1 022 IF (K . EQ. 1 ) GO TO 1 024 _ 

RTHZ (K) = RTHZ (K-1) 

GO TC 102 5 

10.24_STHZ (K) = 1. 

1025 TEMP=(-NBZ*LBRZ*TIMEZ (K)-2.*LWZ*TIMEZ (K) ) *1.0E-06 
RTH1Z (K) =RTHZ (K) *EXP(TEMP) 

RTH2Z (K) = EXP ( ( (TEMP/1 . QE-0 6) -NSCZ*LSCZ*TI M/3600. ) *1. OE-06) : 

C PROPELLANT TANK BLADDER RELIABILITY 

1023 IF(MODEl.EQ.2) GO TO 1026 

RB.1Z (K) =EXP(-NPZ*LBZ*TIMEZ (K)*1 .0E-06) 

RE2Z (K) -RB1Z (K) 

GO TC 1027 ! 

10 26 I F ( MCDE2 . EQ. 2) GO TO. _1 028„ •« 

RB1Z (K) = (EXP (-NPZ*LBZ*TIMEZ (K) *1. 0E-06) ) *ERFC21 

RB2Z (K) =EXP (-N PZ*LBZ*TIM EZ (K)*1. 0E-06) 

: GO . TO 1 02 7 ; • . ■ • . 

1028 RB1Z (K) = (EXP (-NPZ*LBZ*TIHEZ (K)*1 .0E-06) )*ERFC21 
RB2Z (K)=RB1Z (K) 

— C- .RELIABILITY OF PRESSURANT AMD PROPELLANT TANK SOLENOID VALVE SYSTEM 

1027 IF(VSP. EQ. KS) GO TO 1029 
IF (VSP. EQ.KDS) GO TO 1030 

IF(VSP. EQ. KDP) GO TO 1031 

IF (VSP. EQ. KQ) GO TO 103 2 . ~~~ 

C QUAD CONNECTED VALVE SYSTEM 

RSP1Z (K) = (1 . - (1 . -RC1 Z) * (1.-RC1Z) ) **2- ( 1 . - RO 1 Z*RO 1 Z) **2 

RSP2Z (K)= (1.- (1.-RC2Z) * (1.-RC2Z) ) **2- (1 . -R02Z*R02Z) **2 
GO TO 1033 

__C SINGLE VALVE SYSTEM L \ 

1029 RSPIZ(K) = R01Z+RC1Z-1. _ ~~ r 

RSP2Z (K) =R02Z +EC2Z-1 . 

G O TO 1033 . 
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C DUEL SERIES VALVE SYSTEM 

1 030 RSP1Z (K) =RCl.Z*RClZ- (1 .-ROlZ) * (1 .-ROlZ). 1 

RSP2Z (K)=RC2Z*RC2Z- (1 .-R02Z) * (1.-R02Z) 

GO TO 1033 

_£_D.U E L_ P. A R A L L E L VALVE SYSTEM 

1031 RSP1Z (K)=R01Z*R01Z-(1 .-RClZ)* (1.-RC1Z) 

RSP2Z (K) =R0 2Z*R0 2Z- (1.-RC2Z) * (1.-RC2Z) 

GO TO _ 1 0 3 3 : 1 

C QUAD VALUE SYSTEM 

1032 RSP1 Z (K)= (1.- (1 .-ROlZ) * (1.-RC1Z) ) **2- (1 . - BC1Z*RC1Z) **2 

RSP.2Z (K) = (1 . -( 1 .-R02Z) * (1 .-R0 2Z) ) **2- (1 .-RC2Z*RC2Z) **2 — 

C RELIABILITY OF WELDED CONNECTIONS BETWEEN PROPELLANT TANK 
C- SOLENOID VALUE AND THRUSTER 

10 33 RLTZ (K) =EXP (-N WTZ*L WZ*TIMEZ (K) *1. 0E-06) : : 

C RELIABILITY OF WELDED CONNECTIONS BETWEEN PROPELLANT TANK 

.C AND SOLENOID VALVE SYSTEM 

______ R L P P Z ( K ) = E X P ( - N W P Z * L W Z *T . I M E Z { K )*1.Q E - 0 6 ) 1 

C GENERAL SYSTEM RELIABILITY 
1003 IF (SPFVZ. EQ. KN) GO TO 1034 

R FV 1Z=1 . 

RFV2Z=RFVZ (K) 

GO TO 1 035 

10 34 R FV 1 Z_= R FV Z (K )**NPZ ; 

RFV2Z=1. 

1035 IF (CCNVZ. NE.KSQ) GO TO 1036 

B Z_(K ) =R S V 2 Z ■ T 

GO TO 1037 

1036 BZ(K)=RSOLCZ 

_CL_R E PR ES SU R I ZA TI ON NOT REQUIRED FOR EITHER HALF-SYSTEM __i 

1037 RE1Z (K) = 1.- (1. -RFV2Z* RPD 1 Z (K) *RF1Z (K) *RST1Z (K) *RTH 1 Z ( K) * 

1 RLTZ (K) ) * (1 .-BZ (K) *RFV2Z*RFD2Z (K) #RF2Z (K) *RST2Z(K) * 

2RTH2Z (K) *RL.TZ (K) ) . . _ . . 

RE2Z (K) =1 (1. -BZ (K) *RFV 2Z*RFD1Z (K) * R F 1 Z ( K ) #RST1Z (K) *RTH1Z (K) 

1*RLTZ (K) ) * (1.-RFV2Z*RPD2Z(K) *RF2Z (K)*RST2Z (K) *RTH2Z (K) * 

2 RL.TZ (K) ) . : : 

C REPRESSURIZATION REQUIRED FOR FIRST HALF-SYSTEM 
R1Z=RB1Z (K) *RFV1Z*RSP1Z (K) *RLPPZ (K) 

__R2Z = RB2Z (K) *RFV2Z*RSP2Z (K) *RLP.PZ (K) * 

C REPRESSURIZATION REQUIRED FOR BOTH HALF-SYSTEMS 

Q 1 Z= 1 . -fil Z 

Q.2Z = 1 . - B 2 z ; : : 

IF (MODE1 . EQ.2) GO TO 1038 
DEL 1Z= 1 . 

DEL2Z = .1 . • ; 

GO TO 1039 

1038 I F(MCDE2. EQ. 2) GO TO 1040 

DELI Z = R FVZ (K) L*RTZ (K) *RS V 1 Z : 

DEL 2Z = 1 . 

GO TO 1039 

1 040 DEL 1Z =R FVZ (K) *BTZ (K) *RSV1Z _ ■ 

DEL 2Z=DEL 1 Z 

C PROBABILITY OF HAVING DEFFERENT LEVELS OF PROPELLANT 
_C AVAILABLE WHEN NEEDED AND FUNCTIONING PROPERLY 

1039 GZ = R LHP Z (K) *RF V Z (K) *R PTZ (K) 

PFZ (K, 1) = ( 3 . * G Z ) * (R1Z*Q1Z*Q’1Z*RE1Z (K) *DEL1Z* (1.-GZ* 
1 DEL2Z+GZ*DEL2Z*C2Z*02Z*02Z) ♦ B2Z*02Z*RE2Z m * ' _ 
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2 (1. -*GZ *DEL 1Z+GZ*DEL1Z*Q1Z*Q1 2*Q1 Z) ) 

PPZ (K,2)= (3. *GZ) * (R 1 Z*R 1Z*Q1Z*RE1Z (K) *DEL 1Z* ( 1.-GZ*DEL2Zt_ 

1GZ*D£L2Z*Q2Z'*Q2Z*Q2Z) +R2Z*R2Z*Q2Z*R E2 Z ( K) *DEL2Z* 
2(1.”GZ*DELlZ+GZ*DELlZ*QlZ*QlZ*QlZ)+3. *R 1Z *R2Z*Q 1 Z 

3*LQ1Z*Q2Z*Q2Z*RE1Z (K) *RE2Z (K) *DELlZ* CEL2Z*GZ) 

PFZ(K, 3) = GZ* (RlZ*RlZ*RlZ*RElZ (K) *DEL1Z* { 1 .-GZ*DEL2Z + 
1GZ*DEL2Z*Q2Z*Q2Z*Q2Z) +E2Z*R2Z*R2Z*RE2Z (K) *DEL2Z* (1.-GZ 

2*DEL1Z+GZ*D£L1Z*Q1Z*Q1Z*Q1Z) +9.*R1Z*R2Z*Q1Z*C2Z 

3*RE1 Z (K) *RE2Z (K) *DELl Z*DEL2Z*GZ* (R1 Z*Q2Z+Q 1 Z*R2Z) ) 

PFZ (K, 4) = (3.*RE1Z (K)*RE2Z (K) *GZ*GZ*D EL 1 Z* EEL2 Z*R1 Z*R2Z) 

1* (R1Z* R1Z*Q2Z*Q2Z + C1Z*Q1Z*R2Z*R2Z +3.* R 1 Z * S 2 Z*. ; , 

2Q 1Z* Q2Z) 

PFZ (K,5) = (3 . *R El Z (K)*RE2Z (K ) *GZ *GZ*DEL 1Z*DEL 2Z*R 1Z*R 1 Z* 

_____ 1R2Z*R2Z) * (R1Z*.Q2Z + Q1Z*R2Z) 1 

PFZ (K,6)=RE1Z (K) *RE2Z (K) *GZ*GZ*DEL1Z*DEL2Z*B1Z*R1Z#R1Z 
1*R2Z*R2Z*R2Z 

_Q PROBABILITY . OF.HAVING . MORE. THAN THE.DESIREE FRCPELLANT_ _____ 

C AVAILABLE WHEN NEEDED AND THE SYSTEM FUNCTIONING PROPERLY 

IF(HALFZ. EQ. 2) GO TO 1041 

PMZ=PM+PMO •. ; 

GO TC 1042 

1041 PMZ=PM 

1042 FSE= (5 . / 6 . ) * P MO* 2 . 

TTE= (2./3.) *PMO*2. 

OTE= (1./3.) *PMO*2. 

OSE= (1./6.)*PMO*2. 

IF (PMZ.GE.FSE) GO TO 1043 

IF ( PMZ. LT. FSE. AND . PMZ . GE . TTE) GO TO 1044 

J.F_(P M Z . L T . T T E . A N D . P M Z . G E . E M 0 ) G 0 T CL._1. 0 4 5. 

IF (PMZ. LT. PMO . AND. PMZ. GE .OT E) GO TO 1046 
IF (PMZ.LT.OTE. AND . PMZ. GE . OSE ) GO TO 1047 

IF (PMZ. LT. OSE. AND.PMZ.GE.O.) GO TO 1 04 8 ; 

WRITE (6, 1049) PMZ 

1049 FORMAT (IX, 'PMZ IS IN ERROR PMZ=',F10.6) 

CALL. .EXIT • : : : 

C AT LEAST ONE PROPELLANT TANK AVAILABLE 

1043 IZI=1 

_______ G 0_.T 0 _ 1 0 5 0 __ ... , 

C AT LEAST TWO PROPELLANT TANKS AVAILABLE 

1 044 IZI=2 , 

GO TO .105 0. • ... - 

C AT LEAST THREE PROPELLANT TANKS AVAILABLE 

1045 IZI=3 

GO TO. 1050 . ; t 

C AT LEAST FOUR PROPELLANT TANKS AVAILABLE 

1046 IZ 1= 4 

GO TC 1050. . : : 

C AT LEAST FIVE PROPELLANT TANKS AVAILABLE 

1047 IZI= 5 

, .GO TC 1050 , 

C ALL SIX PROPELLANT TANKS AVAILABLE 

1048 PZ (K ) =PFZ (K , 6) 

GO . T C_. 1 0 6 5 ' . - 

1050 PZ (K) = 0 . 

DO 1052 IZ=IZI , 6 

_ J.052 PZ (K) = PZ (K) +PFZ (K,IZ) . ... 
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GO TO 1065 

16 ICTR-0 — : : — 

WRITE (6, 54) 

WRITE (6,61) | 

G O TO .,6 : : : — j 

21 IERR=1 . ; 

GO TO 31 

10 65 IF (I ERR . EQ. 1 ) GO . T0...3 1 ; 

40 CONTINUE •, 

3 IF (HALFZ.EQ.2) GO TO 7 

H ALFZ=2 

KI=NK+1 
J= 0 

wsc=.wsc+.pa ' - : 

GO TO 2 
7 I NK= 0 

__1_6_6_.I-C.TR = 0 : * : - 1 

IK=I NK + 1 

WRITE(6,72) ' 

72 FOR MAT (1 HI , 7X, * MANEUVER' , 1 OX, 'TOTAL* , 17X,* BORN* , 20 X, '.DELTA ' ,18X,_ 
1 'TOTAL' ,1 6 X,' TOTAL' ,/,9X, ' NUMBER* ,10X,'It1PUISE',16X,*TIME',22X, 
2'V • , 19X, • DELTA V ' , 13X, • BURN TIME ' ,/, 2 5 X , • (N-S EC) • , 1 5X , ' (SEC. ) ' , 

318X,-'. (M/SEC) • , 1 6 X ,'_ (M/SEC) • , 15X ,2_(S.EC) !,/).„ 

DO 18 K=IK,NK i 

WRITE (6,74)K # STI(K) # STIHE (K) # SDELV (K) ,STDEL(K) ,STIM(K) 

74 FORMAT (IPX ,13, 8X,E12.4, 8X, El 6 . 6 , 8X, E17_. 4 ,8X,E14...4 , 8 X ,F 1 0 . 0 ,/) 

ICTR=ICTR + 2 
INK = K 

I F f IC TR. GS. 54 )_GQ_TCL-6j6 — : 

18 CONTINUE 

C OUTPUT K VARIABLES j 

INK = 0 : , — 

1051 ICTR=0 

IK= INK+ 1 1 

HIT E(6,1053) . : — 

1053 F0RMAT(1H1,4X, 'MANEUVER', 2X, 'TIME', 5X r ’ SYSTEM ', 4X, 

1 ' LEV EL OF SYSTEM AVAILABILITY (NO. OF PROP. TANKS) • ,6X, 

21 RELIABILITY OF ',/, 6 X, ' NUMBER ', 3X ,' (HRS) ' , 4X, ' RELI AB ' , ; — : 

36X,'1*,8X,'2',8X,'3',8X,'4',8X,'5',8X,'6' ,7X, 'THRUST SYSTEM' , 

4/, 8 X , ' K' , 6 X , ' TIMEZ' , 6 X , ' PZ • , 4X , • PFZ (K , 1) ' , IX, 'PFZ (K,2) ' , 

5 IX,. 'PFZ (K ,3) ', IX, 'PFZ (K, 4) ', IX, 'PFZ (K,5) • , 1 X , ' PFZ (K , 6 ) ».,_4X, 

6'RE1Z (K) • ,1X,'SE2Z (K) ' ,/) 

DO 1054 K=IK , NK 

WRITE (6,1055) K, TIMEZ (K) ,PZ (K) , (PFZ (K, KK) , KK=1,6) , RE1Z (K) , 

1RE2Z (K) 

1055 FORMAT (7 X, IJ,3X,F9. 1 , 7 ( 1 X ,F 8 . 6) , 2X , 2 ( 1 X , F 8 . 6) /) 

ICTR = I.CTR +2 — 

I NK=K 

IF (ICTR .GE.54) GO TO 1051 

1054 C 0 N T I N U E 1 : 

INK = 0 

1056 ICTR = 0 

IK=INK+ 1 ; 1 — 

WRITE (6,1 057) 

1057 F0RMAT(1H1,4X, 'MANEUVER' ,3X, 'NUMBER OF * , 5X , • FILL* ,3 X , 

1_* HI_..PRESS_.PROP_TANK' ,4 X,' PROF TA NK ', 9 X , • BLADDER * , 1 1 X , 
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2' PRESSURE',/, 6X, • NOM BER 5 X, • CYCLES • , 7 X,* VALVES' ,4X, 

3 • I’AN K* , 3 X , ' 51 ELD RING' , 3X, • VALVE SYSTEM ' »5X, '8 ELI ABILITY ' , — 

48X,' TRANSDUCER ’ , / , 8X , ' K ' , 5 X , ' C YCL1 Z CYCL2Z* ,4X,'RFVZ' ,5X , 

5'RTZ* ,7X,» RPTZ* ,5X,'RSP1Z« , 4X , ' RSP2Z ' , 5X, *RB1Z',4X, * R B 2 Z * , 

65X,_'J*ED1Z',.3X, ' RPD2Z ' ,/) : 

DO 1 05 8 K=IK,NK 

WRITE (6, 1059) K , C YCL 1 Z ( K) ,CYCL2Z (K) ,RF VZ (K) , RTZ (K) ,RPTZ (K) , 

IRS PlZ(K) . R SP 2Z (K) ,RB1Z (K) ,R32Z (K) ,RPD1Z (K) ,.RFD2Z (K) . — 

1059 'FORMAT (7X , 13, 3X, 2 ( IX, 16) , 2 F 9 . 6 , 1 X,F 9 . 6 , 1 X , 2F9 .6 , 

12F9.?6/2F9. 6,/) - - - ' 

ICTR=ICT R +2 4— 

I NK=K 

IF (ICTR .GE.54) GO TO 1056 

1058 CONTINUE 1 Li ; 

IN K=0 .. .. 

1060 iCTR-0 * 

■ IK=I N K +1 — — 

WRITE (6,1 061) 

1061 FORMAT (1Hl,4X, 'MANEUVER ',7X,' FILTER • ,10X,' THRUSTER* ,9X, 

1« T HRUSTER' , 6X, » HALF-SYS' ,2X,'HELD CONNECTION RELIAB',/, 

2 6 X , ' NUMBER ' , 7 X , ' ASS EM EL Y ' , 7 X , ' V AL VE SYSTEM' ,6X, 

3' RELIABILITY' ,4X, 'CON. V A'LV 1 , 2X , ' HP SYS', 2X, 'PRO SYS', IX, 

4 'THR SYS',/,8X, 'K',8X,'RF1Z* ,4X,'RF2Z' ,5X,'.RST1Z* ,4X,. 

5* RST2Z' ,4X ,' RTH1Z' ,4X, • RTH 2Z ' ,5X,’BZ» ,6X, ' RLHPZ ' , 4X, 

6'RLPPZ',4X, ' RLTZ ' ,/) 

DO 1062 K=IK,NK ' I 

WRITE (6, 106 3) K,RFlZ(K) ,RF2Z (K) ,RST1Z (K) ,RST2Z (K) , RTHIZ(K) , 

1RTH2Z (K) ,BZ (K) , RLHPZ (K) ,RLPPZ (K) , RLTZ (K) 

JQ6 3_JPOR MAT (7X , I3,3X, 2F9 . 6 , 2F9. 6 , 2F9 . 6 , F9. 6, 

13F9.6,/) 

ICTR=ICTR+2 

I NK=K ; 

IF (ICTR. GE.54) GO TO 1060 

1062 CONTINUE 

CALL_ EXIT 1 

100 WRITE (6,63) 

63 FORMAT ( ' READ IN ERROR') 

C ALL EXIT .. - — — — : r 

114 WRITE(6,68) C3,C1,PG 

68 FORMAT (' FM IS LE 0. , C3 = ' E16 .6, * Cl = ' ,E16. 6, ' PG - ' ,E16. 6) 

CALL. EXIT : * : : : : 

54 FORMAT ( 1 H 1 ,25X) 

END 
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0027 

0028 

0029 

0030 

0031 

0032 

0033 

0034 

0035 

0036 

0037 

0038 

0039 

0040 

0041 

0042 

0043 

0044 

0045 

0046 
004 7 

0048 

0049 

0050 


SOBBCOTINE STINT (A EG1 , A RG2, ABG3 , PC T , KE Y, NGRIPB.NINTBL , K AXTPL) 
DATA K002PX/5 / 

DATA K001FK/6 / 

DIHENSI OSNUN PTS (38) ,L1 (37) ,L2 (37) ,L3 (37) ,STG (1765) ,DUNNY (10) 
DIHENSION NAME (9) 

C DINE BSICN NO NETS ( 38) , Li ( 3 7) , L2 (37) ,L3 (37) , STG (1765) ,DOSN Y ( 10 ) 

C DIHENSION NAHE (9) 

EQUIVALENCE (NAT,L3(1)) 

C***« SIZE OP STG IS CALCOLATED BY SON OP ( (1 ♦» (A BG1) ) * ( 1 »N ( APG2) ) ) 
NGHIPE=0 
I F (KEY) 1,1,70 
1 NG= 1 
N0RNAL=1 

WRITE (KO0 1? X, 1 357 ) 

1357 FOB a AT (34H0 TABLE CATE CONTENTS) 

GO TO 55 
2000 NG=2 

MORHAL— 2 
3000 RETORN 

775 NGRI EE= 1 
RETOEN 

776 NG8IPE=2 

WRITE (KOO IPX, 9000) AHG1 , ARG2 , ARG3 , SINTB l, KAXTBL 
RETURN 

9000 PORNAT (20HO ERROR IN TLO , ARG 1=F 1 2. 5, 6H A RG2=P1 2 . 5 ,6H ARG3 = F12.S 
18H NINTBL=I4, 8H HAXT3L=I4) 

C GRUNNEN AIRCRAFT ROUTINE FOLLOWS 

C BEGINNING OF STINT 

• 55 NUHTEL=1 

NONPTS(I) =0 

102 -READ (K002FX, 57) DAI, CA2, DA3,K,L1 (NOHTBL) ,L2 (N0H7BL) , SAHE.ISEQ 
57 FORHAT (A2.A3, A 3. 1 4, 21 2, 9 A 4 , 18X.I2) 

WRITE (K001FX.1157) K , CA 1 , E A2 , C A3 , N AHE 
1157 PORHAT (IR,5X.A2,A3.A3,5X,9A4) 

104 IP(ISEQ) 69,58.69 

59 I? (K) 99,99,1159 

1159 IF (K-37) 59, 59,1103 

59 LS=L1 (NUHT3L) 

N 1= (L8- 1 ) /9 ♦ 1 
DO 60 IS — 1 , N 1 
NAT= (IS-1) «9*1 

IF (IS-N1) 60,61,60 

60 L4=NAT*3 
GO IC 62 

61 L4=L8 

62 L5=NU.NP?S (NUHTBL) *1 
L6=L5>NAT 
L7=L5»L4 

JJ=0 

L9-L2 (NOHTBL) 

LN=L5+L3 

LN=LH>L9 

105 READ (K002FX.64) (DUHHY(K) ,K=1,10) ,ISEQ 
WRITE (K001FX ,44) (DUNNY (K) ,K=1 ,10) ,TSSQ 

44 PORHAT(10S12.4,3X, 12) 

64 PORNAT (1CE7.0, 12) 


01 C0010 
01 00020 
01 C0C30 
01 00040 
01 0005C 
01 00060 
01 C007C 
01 00080 
01 00090 
01 00100 
01 00110 
01 C012C 
01 00130 
01 C0 14 0 
01 00150 
01 C016C 
01 00170 
01 001R0 


01 C0230 
01 00240 
01 00250 
01 C0260 
, 01 00270 
01 CC29C 
01 00290 
01 00300 
01 00310 
01 CC32C 
01 00330 
01 C0340 
01 C0350 
01 C0360 
01 CO 370 
01 C0380 
01 C0390 
01 00400 
01 CC410 
01 00420 
01 C0430 
01 00440 
31 C0450 
31 C0460 
01 C0470 
01 C0480 
01 00490 
01 C0500 
01 00510 
01 00320 
01 C0530 
01 00540 
01 


01 00560 
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0051 

107 

STG(L5)=D0HSY(1) 

01 

00570 

0052 


K=2 

01 

C0580 

0053 


00 65 J=L6,L7 

01 

00590 

0054 


STG(J) =D0MMY (K) 

01 

C0600 

0055 

65 

K=K»1 

01 

00610 

0056 


IF (ISEQ- ((IS-1)*(L9»1)»JJ*1)) 69,66,69 

01 

00620 

0057 

66 

L6=LN*NAT 

01 

0063C 

0058 


L7=LS*L4 

01 

00640 

0059 


L5=LH*1*JJ 

01 

C0650 

0060 


IF (JJ-L9) 67,68,69 

01 

00660 

0061 

67 

JJ*JJ«1 

01 

C0670 

0062 


LN=LN*L8 

01 

C0680 

0063 


GO TC 105 

01 

00690 

0064 

68 

CON TIN OE 

01 

C0700 

0065 

109 

LEB=SOMPTS (NOMTBL) *(L8*1) * (L9*1) 

01 

00710 

0066 


IF (LEE-1765) 1100,1900,1101 

01 

00720 

0067 

1100 

IF (NOMTBL-37) 1102,108,1103 

01 

00730 

0068 

1102 

N0HPTS(NOSTBL* 1)=LE£ 

01 

0074 0 

0069 

108 

SO HTEL=N0 MTBL*1 

01 

C0750 

0070 


GO TO 102 

01 

00760 

0071 

1101 

SBITE (K001FX,1111) LEE 

01 

C0770 

0072 


GO TO 775 

01 

00760 

0073 

1103 

SHITE (K001PX, 1 113) NO NT 8 L 

01 

C0790 

0074 


GO TC 775 

01 

00800 

0075 

1111 

FORMAT (17H TOO MANY POINTS 18) 

01 

00810 

0076 

1113 

FORMAT (1 7H TOO MANY TABLES 18) 

01 

C0820 

0077 

69 

GO TO (775,776,776) ,SG 

01 

00830 

0078 

70 

IF (HINTBL-MAITBL) 71,100,69 

01 

C0340 

0079 

71 

DO 73 NAT=MINTBL, MATTEL 

01 

00850 

0080 


L4=NDMPTS (NAT) *1 

01 

00860 

0081 


IF (ARG3-STG (L«) ) 72,74,73 

01 

00870 

0082 

72 

IF(NAT-MINTBL) 69,69,75 

01 

C0880 

0083 

73 

CONTINOE 

01 

C0890 

0084 


GO TO 69 

01 

00900 

0085 

75 

L5=1 

01 

C0910 

0086 


L6=2 

01 

00920 

0087 


L7=L4 

01 

C0930 

0088 

101 

DO 97 L8=L5,L6 

01 

C09 40 

0089 


L4=H0HPTS (NAT) ♦ 1 

. - 01 

00950 

0090 


L9=L1 (NAT) 

01 

C096C 

0091 


LM=L9*L4 

01 

00970 

0092 


DO 77 LN=1,L9 

01 

00980 

0093 


JJ=L4>LN 

01 

00990 

0094 

2626 

IF (ARG1-STG (JJ) ) 76,78,77 

01 

01000 

0095 

76 

IF (LN-1) 69,69,79 

01 

01010 

0096 

77 

CONTINOE 

01 

01020 

0097 


GO TO 69 

01 

C 1030 

0098 

78 

N1-- 1 

01 

01040 

0099 


GO TO 80 

01 

01050 

0100 

79 

N 1 = * 1 

01 

01060 

0101 

80 

K=L2 (NAT) 

01 

01070 

0102 


DO 82 1=1, K 

01 

C1080 

0103 


IDATE=LH*I 

01 

01090 

0104 


IF (ARG2-STG (IDATE) ) 81,83,82 

01 

C1100 

0105 

81 

IF (1-1) 69,69,84 

01 

dnio 

0106 

82 

CONTINOE 

01 

C1120 
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0107 

GO TO 69 

01 

C1130 

0108 

83 IS=- 1 

01 

01140 

0109 

GO TO 85 

01 

01150 

0110 

84 IS=*1 

01 

01160 

0111 

85 ISEQ=LS*L2 (NAT) ♦LN* (1-1) *L9 

01 

Cl 170 

0112 

J=IS£Q-L9 

01 

01180 

0113 

K8=LH* (1-1) 

01 

01190 

0114 

K9 a I4*LN— 1 

01 

Cl 200 

0115 

, IP (Nl»IS) 66,88,91 

01 

01210 

0116 

86 IP (STG (ISEQ) -999. E20) 87,69,69 

01 

01220 

0117 

87 FCT=STG (ISEQ) 

01 

01230 

0118 

GO TO 95 

01 

01240 

0119 

88 IP (HI) 89,69,93 

01 

01250 

0120 

89 IP (ANAX1 (STG (ISEQ) ,STG(J))-999.E20) 90,69,69 

01 

C1260 

0121 

90 FCT=STG (ISEQ) -(STG (I DATE) -ABG2) * (STG (I SEQ) -STG (J) ) / (STG (ID ATE) 

01 

01270 


1-STG (K8) ) 

01 

01280 

0122 

GO TC 95 

01 

C1290 

0123 

91 IP (AHAX1 (STG(ISEQ) , STG (J) ,STG(ISEQ-1) ,STG (J-1 ) ) -999. E20) 92, 

01 

01300 


. 169,69 . 

01 

01310 

0124 

92 PCT= ( (STG (ID ATE) -AHG2) * ( (STG (JJ) -ABG1) *STG (J-1 ) - (STG (K9) -AEGl) 

01 

01320 


1*STG (J) ) - (STG (K8)-ARG2) * ( (STG (JJ)-ABGI) »STG (ISEQ- 1) -( STG (K9) - 

01 

01330 


2ABG1) *STG (ISEQ) ) ) / ( ( STG ( IE ATE) -STG(K8) ) * (STG (JJ) -STG (K9) ) ) 

01 

01340 

0125 

GO TO 95 

01 

Cl 350 

0126 

93 IF (AHAX1 (STG(ISEQ) , STG (ISEQ-1) ) -999. E20) 94,69,69 

01 

01360 

0127 

94 PCT=STG ( ISEQ) - ( STG ( J J) - ABG 1 ) * (STG (ISEQ) -STG (ISEQ-1) )/ (STG (JJ) - 

01 

Cl 370 


1STG (K9) ) 

01 

Cl 38 0 

0128 

95 50 TC (96, 98,99) ,L8 

01 

01390 

0129 

96 DUNS T (1 ) =PCT 

01 

Cl 400 

0130 

97 NAT=HAT-1 

01 

01410 

0131 

98 PCT=DUnBT (1)- (STG (L7) -ARG3) * (DUBBY ( 1) — PCT) / (STG (L7) —STG (L4) ) 

01 

Cl 420 

0132 

99 GO TO (2000, 3000) .NOBBAL 

01 

01430 

0133 

100 HAT=BIHTBL 

01 

01440 

0134 

74 L5=3 

01 

01450 

0135 

L6*3 

01 

01460 

0136 

GO TO 101 

01 

C1470 


C EHD STINT TABLE LOOK-UP 

01 

01480 

0137 

END 

01 

Cl 490 



B-l 


APPENDIX B: 


Printout of Input Data and Computed 


Results 


^ m 

iurff a 
0000 
y m y y 
• 0 © O' if 


• 0*00000000 
■w 

OOOOOOOO' 


O O O' 

U 0 0 
OW« 

• 0 X 

© 0 • 

O 0 O 0 

*«» 0 r- 0 


f *n VI 
000 

0 0 0 
*n 0 X 
O' p * it 


• to 


<d oooooooo 

aoooooooo 


B O 9> O ffl O' 

00000 

O 0 O' ^ *0 
• O' 0 0 0 0 
O O' 


».«- o o' 


• *-• • X 
OHS' 
• at O' 

M O* 


• o O' 

> . O* 

o o « 

• *<N 

«o r* c 
r» r- < 
o *cn c 


• o'r-ooooooof- 

• * 

© O O O O © O 


» O' » 
> O' O' 
' O' O' 
'O' O' 
' O' » 
4 04 tO 
• ® 

■ • 

3 r- r» 


I** a 
O' O' 
0 o> 


xr-ooooooo 

• • 

o o © © o o © < 


o o 1 

• • o < 

« o o 
m o r 
* *o < 


‘ O* O' 
O' o* 
O' O' 
\ o* o* 

■ a p- 
i o « 


tfOS « 

O' o O' 

O' if O' 

o* f-‘ \0 . C 

O' o o* 

O' • O' 

«n «0 • 

p-r-oooo©oor 

p- m . • 

o o o o o o o 


x.0 

O' O' 

I IN O' 

) 03 O' 


o o 
• • • 
o o 


» CO O' : 

> O' O* c 

> O' O' c 

> *f 0 *. 

>»■'«( 

> *o a • 

» • • 

• tn f 


r-ooooooo • 

-- 

oooooooo 


M*o m 0 0 0 0 r- 


oo i 
I • ♦»* 
a © o a, i 


x m « 

O* O' c 
O* O' c 
0 0 C 
oo r 
o< in i 
in r* 


• 0*0000000 

“'*••••••• 

o o o o o o o 


,fci; 


' O' O' 
' O' O' 
' O' O' 
' X O' 


• o o 
© o 


O' O' 
O' O' 
O' O' 
O' O' 

»n o 
o r- 


> •- » 

• o • • • 

• 0 p* ic n 


<X> O' 

O' ® 

O' O' 

O* O' 

o* *- 

»n in 

•XOOOOOOO 

n in • 

o o o o o o o 


0000 m ( I 

0*000 riN N 

crooo x (N ca tt 

if O O O •CO 4N 

oooo p- x n » » 

• • • • • O «N X 

oo*» * i »• a. a 

or**xino©ooo©ow 0 .-i ^ 


o o o o o o 


> o o o* 
» o O if 
>003 


O O' ui o 
0*0 O 

*o* 0 o 

o* m o 

JnC 003 

O' x o o o 

31/1000 

o r- o *- o 


m t m * o 


u » **» t 
•4 N C 

>0 % > J 


0030000 


> O O O 
>030 
3 0 0 0 
3 0 0 0 

> O o o 


O' o 
O' o o 

o *n 

• O N 

mo • 
o o 


0 * o 

O' O 

O' o 


» O • O O' o 
• • • • m *- r- u • 
>ooo» m cn o 


• o o> o O' 
s o c a r* 
if o O' o o* 

© if © if 

o o o o 

• o • • • 

irrNfCNQOOOOOc 




OOOOOOOOP 


in © m 

o* © rvt 

o* © o • 
O' .00 

7*fVOO 


> o o o 

>03© 

> O O © 

>03© 

> o o © 

> • • • 

• O' 0 at 

r'NrsifNOO©oo< 
; cm f • • • • • 
o o o o o c 


• O' o 


• if © 0 0 
f y o o* p* 
m © o © © 

a> o nm 
f o o o 

• o • • • 

• <e o 

ON «-^000000 
pr- mm 

o o o o o o 


w © m 

*- ©in 

'•KOI 

• o M o 

> X © * 

■ X O 0 *. 
IN * O l/i 
■a in w~ o » 
r> » © O • 
*-•-> I o 
l »oa • 

• *■ © o 

- o © «- 

I © o 
© in 
© © « <n 
O' O *4 • 
© O tJ 0 
9«NPO 

o* • cn f- 

0 © »in 



f <N 4f 0 

**NU» 

pNff » 

P- IN 













r* ■ 






M 

O <N 

r* 0 

r- *“ 

IN 3 


X2 

© © 

o o 

3 o 

9 O - 



Ul Ik] 

Ul *1 

Ul M 

M X 



© O 

IN X 


o O 



o o 

X if 

if 0 

o © 



r- f 

0 N 

0 sj 

IN 30 


•f 

0X00 

•-moo 

•-*-0 0 

«- 0 o o 

© o 




• • • • 

■ • « • 



oooo 

3 0 0 0 

oooo 

oooo 

© o 

■4* 






X 

O fN 

r- 0 

r* — 

N © 



o o 

© © 

o o 

3 3 




ul 1) 

Ul w 

ui X 



o o 

0 03 

0 0 

o o 



o o 

fn -N 

IN 3 

o r— 


X 

*o o 

•f m 

r- P* 

o o 


-4 

0 in o o 

•- CN O O 

•- — © o 

*- 0 O o 






• • • • 


» 

oooo 

OOOO 

oooo 

OOOO 








cc 






o 

O IN 

r* 0 

r- p- 

p* o 






o o 


X 






O 

(U X 


CJ M 

M N 


•4 

o o 

»- X 

if o* 

© o : 


u. 

o o 

0 O 

x *- 

O o 


HI 

m m 

IN IN 


O \C 


a 

0 N O O 

— IN O O 

4- P- © O 

X 0 © o 

o o 

O 






l 

oooo 

oooo 

oooo 

O ooi o 

o o 

m 






X 

O p- 

3 0 

r- •- 



■a 

o o 

O o 

o o 

o © 


o 






u 

w rq 

M M 

03 Ul 

a ui 1 



o o 

0 ao 

if in 

o o 


X 

o O 

•N ® 

r- in 

3 © 


u 

0 © 


*- r* 

O 0 


HI 

0 0 O © 

*- m O O 

•- — © © 

*o r- o o 

o o 

a. 




• * • • 

• • 


oooo 

oooo 

OOOO 

o o o o 

3 O 

CB 






O - 

a i 

© *- 

r* 0 

■0 ^ 

p- © 


a t 

© © 

o © 

o o 

© © 


u ; 

OJ X 

M Ul 

03 U 

M X 


« i 

o © 


-N *43 

© © 


Ml 

o o 


IN — 

© © 


to 1 

f in 

a p- 

X *D 

O 0 


X ' 

0u>oo 

•-moo 

X — O O 

m r*. o © 

3 © 

3 1 


• • • • 

• • • • 

• • • • 


X 

© o o © 

© o o © 

oooo 

© o o © 

© o 

X . 




p- o 

if if 

S! 

o o 

© o 

. o o 

3 © 

© o 

1 


M 


m 


QJ X M 

E © O 

0.0 0 

f rvi © 

M 0 © O o 

u • • • • 

© o o © 


BOO 


O o • 
• • IN 

o o o 


• N 

o u * 

IflOH 
n o in 
m •© » 

Cu o • 


x a a 
m © 
rN f 

«c «- 

cd rN o o 

©003 

HI X 0 © 0 ► 
H4 © O O © l 

J U Hi N W . 
PP p OM 

* y B Oin : 

M X © © CN * 

0 in 3 © ifl (N p N 
U • • • • z • • • •; 

0.000003000 : 

r 

w 

» O p- ^ rr I 

O © © © i 


ttl 4) 

r- 0 
© m 
© m . 

sC p- O o 
• • • • 
© Oo o 


0 X © O » 

• • • • c 
oooo; 


3 u x 0 r- 0 ( 
> a. o © © o : 


» O' 
X © 
*-* 0 
K 0 


• oooo 

p- O C-- O o 
fM © O © © 

©003 

© O © O 

• 0 0 

©•“•-•-OOOOOO 
P p fNPI ..... . 

O O O O O o 


»MPOO 

©oooo 
0 i o o m 

O tu O © IN 

0 0 o o • 
0 0 o o o 

0 0 O O <N 


. Hi Z M u j] 
s Q © © © 3 
4 OOOO 
* © © O 3 

s 0 m ^ in 


WhiUUWl 
» T IN 0 <N : 
rn 0 cd in 
PJ *- 0 3 IN ( 
W 3 p p IN ! 


OOOO oooo 


W Ml 
m f 
P» *45 

r- 0 

in •- © © 

oooo 


u z 
o © 
r» in 

0 P- 

I N pp © o 
• • • • 
oooo 


JC OJ M w tt* 
3 0 0 0 0 
X O IN © O 
0X30 
Mnmpp 
U5 


f o *n 
on.© 
O m O' 

a in «- 


uooooao© 


U U !il ill 

B 3 O © O 

tw © J O o 

w o in © © 

IN 0 N P 


M U 
O 0 
© O 
IN O 


OOOO o o 


O O II 
• . N 
O O X 


It i£> 

in in < 

cu p- : 


M 1-4 
OttO 


VOUU HOirvt4«V4* 


o ••••••••*• 

© 

© 

o 

o 

• oooo 

o © © © © 

<n o © © © 

o o © o 

5 © O © © 

■« II O • • • 

e N • f 0 © 

~ M0«- 00000300 

C p - p- f • • • • • • 

HI © © © O © © 


P- II 

• n n n 

O CJ N (J 
in x m 
ii z r > 


• N 


© © © OOOO 
p- UWWpMMWUp 

r- ©o©r**o©vD®f»* 

W | 0 0 3 | (N *0 ® 0 | 

Hf ©Ein0r‘Xtn<N0 

<OOmc»* l NpppNNO 


UZUpZZ 

3 © *- r- o o 
3 O I O O 
0 0 0X30 


o I 


H-OOOOINOOOOINOOOOINOO 


xC r— 

p- 0 

© o 


r* o \c 

I N IN 

m 0 o 

nap 


i • • 

moo 


o © 
© © 


© ^ m © «f 
e r on 
© u*o® 
o 0 o o 

O 0 © CD 

© 0 o o 

O 0 o o 

• IN • If O 


woo IA.UI 

4irm«nu 


OOOO 
• • • • 
oooo 


© © © © 
oooo 


© © © © 
oooo 


oooo 
• • * * 
oooo 


o © 

© o 


ViW ilWNl 
BMWdOiBM 


O N 

— 

<J M IN 


N — 
u a h 

* H w 

M » r 

3 

OJ Ml u 

* to — 

M -i 

*• » p» 

a* o 

n# a X 


f o a. in 


M *• *4> 
HUH 


a. in © in 


\« ’ 

H E • 
a m i 
B *- 


o m 
o o. 

•» 

f- i| ft. H 


3 -» — U 

y> O m cn 

-H = B * 
t- 04 u w 
w e \a c 
x won 
3 HP ax u 
© cj a. f 
w w in n -4 
\ x a a* 


LIST or VARIABLES WITH CORRESPONDIN'? UNITS 



B-2 



<M 
' O 
1 

u 

M 

B H ft) 

OV 

ohm 

m 



ft. 0) O 

r* 

kS 

• 


o 




00 



• 

cu 

u 

<T* 

W 

M 

Psl 

M 

tO 

CM 




9 

o 

w 

PM 



<M 


U 

in 


M 

CM 

M 

V) 

CM 

> 

V 

• 


S3 

O 



< 




r- . 

o 





o 

t 



H 






W 

K 

M 

ft 



■1 

r to 

r- 




*< M Ot 

CM 

o 



o» 

S 03 tO 

r* 

r* 



to 

Uft \ 

T* 

ft 



\ 

X 

• 

0- 


r* 

SB 

w 

o 

>> 



r* 




4 


o 






M 






CM 






9 


r* 




CM 


o 






s 

• 

B- 

M 



o 

E CO 

in 


o 


25 M O i 

00 

w 

to 


-< OS to 

in 

E 



H 0* V 

m 


m 

II 

9E 

• 

o 

o 


w 

o 



CO 

1 

z 





M 

o 




MP 

rn 

M 




O 

in 





1 

o 

«c 




fltt 

in 

64 




ar 

• 

W 


m 


o 

o 

03 


o 


o 


S 




CM 


to 


M 

CO 

• 

n 

CO 

CU 6-1 #-» 

r* 

B 

o 


w 

O 0u O 

in 

o 


03 

03 

03 M Z 


M 


u 

04 

CU ftj w 


Eh 

If 

!H 

W 


• 

W 


BLl 

03 


o 

O 

«< 

E 




sb 

M 


E-* 



o 

03 

M 

«< 



u 

< 

Q 







W 

03 


*1 

W 

W 

03 

W 03 



u 

u 

D 

> W 

**- 

H 

w 

M 

CO 

3 O 



04 

fe 

to 

u e 


w 

M 

M 

H 

z a 


SB 

OS 

03 

03 

< z 


M 

o 

O 

C4 

e 




B-3 
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CHltt 
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(H/SQ H) 

THEOST 

(D 


TBJ 
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ISP 
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FI04 
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(KG/SEC) 

1 
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03 

0.33398 
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0.16628 

07 

0.15908 
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0.225108 

04 
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1 * 

0.35418 
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0.76718 

06 

0.56188 

06 

0.46698 

01 

0.203961 

04 


208.13 

0.22898-02 

15 

0. 33558 

02 

0.75308 

06 

0.55331 

06 

0.45848 

01 

0.203698 

04 


207.85 

C.2250E-C2 

16 

0.31988 

02 

0.74158 

06 

0.54638 

06 

0.45148 

01 

0.203478 

C4 


207.62 

C.2218E-02 

17 

0.30048 

02 

0.72778 

06 

0. 53798 

06 

0.44308 

01 

0.203208 

04 


207. 34 

0.21808-02 

18 

0.27728 

02 

0.71198 

06 

0.52818 

06 

0.43338 

01 

0. 202888 

04 


207. C2 

0.21368-02 

19 

0.85148 

01 

0.60308 

06 

0.45588 

06 

0.35438 

01 

0.200568 

04 


204.65 

0. 17678-02 

20 

0.65858 

01 

0.59388 

06 

0.41598 

06 

0.31268 

01 

0.199288 

04 


203.35 

0. 1569E-02 

21 

J). 625.18 

01 

0.59228 

06 

0.39318 06 

0.28988 

01 

.0, 198098 

04 


202.13 

C, 1 462E-C2 

22 

0.11518 

03 

0.34898 

07 

0.17148 

07 

0.16378 

02 

0.225228 

04 


229.81 

0.72698-02 

23 

0.94038 

02 

0.18248 

07 

0. 10768 

07 

0. 10178 

02 

0.217728 

04 


222.16 

0.46708-02 

24 

0.93388 

02 

0. 17978 

07 

0.10658 

07 

C. 10048 

02 

0.217468 

04 


221.90 

0.46188-02 

25 

0.73528 

02 

0.1236B 

07 

0.81568 

06 

0.72938 

01 

0.211658 

04 


215.97 

0.3446E-02 

26 

0.72908 

02 

0.12248 

07 

0.80978 

06 

0.72318 

01 

0.211488 

04 


215.60 

0.3419E-02 

27 

0.71388 

02 

0.11958 

07 

0.79521 

06 

0.70818 

01 

0.211C7E 

04 


215.38 

0. 3355E-02 

28 

0.70858 

02 

0.11868 

07 

0.79048 

06 

0.70318 

01 

0.210938 

04 


215.24 

C. 33331-02 

29 

0.69988 

02 

0.11708 

07 

0.78258 

06 

0.69498 

01 

0.210718 

04 


215.01 

0.32988-02 

30 

0.6S1SP 

02 

0. 11398 

07 

0.76648 

06 

0.67848 

01 

0.21025E 

04 
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0.32278-02 

31 

0.64748 

02 

0.10858 

07 

0.73808 

06 

0.64948 

01 

0.209448 

04 


213.72 

0.31 008-02 

32 

0.36018 

02 

0.77188 

06 

0.56468 

06 

0.46978 

01 

0.204C56 

04 


209.22 

0.23028-02 

33 

0.29688 

02 

0.7253E 

06 

0.53648 

06 

0.44158 

01 

0.203158 

04 


207. 29 

0. 2 173E-02 

34 

0.21828 

02 

0.67468 

06 

0.5055B 

06 

0. 40928 

01 

0.202168 

04 


206.28 

0.20248-02 

35 

0. 15098 

02 

0.63648 

06 

0.48238 

06 

0.38318 

01 

0.201418 

04 


205.52 

0. 19028-02 

36 

0.80958 

01 

0.60108 

06 

0.45 37 E 

06 

0.35218 

01 

0.200498 

04 


204.59 

0.1756E-02 

37 

0.62518 

01 

0. 5922B 

06 

0.39318 

06 

0.28988 

01 

0.19809E 

04 


202.13 

0. 1462E-02 



B-4 


BAB BOV ER 

■ORBEB 

total 

IMPULSE 
(»-S EC) 
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TI8E 

(SEC.) 


DBLTA 

? 

(R/5EC) 

TOTAL 
DELTA V 
(M/SEC) 

TOTAL 
BOBS TIME 
(SEC) 

1 

0.3808B 04 

0. 234000E 03 


0.7006E 01 

0.7006E 01 

234, 

2 

0. 81431 04 

0.2800C0E 03 


0.7807E 01 

0. 148 IE 02 

514. 

J 

0.1288S OS 

0.323000E 03 


0. 8633E 01 

0.2344E 02 

837. 

4 

0.5098B 05 

0. 3239C0E 04 


0.7085B 02 

0.943CE 02 

4076. 

5 

0.5S10E 05 

0.422CC0E 03 


0.7855B 01 

0. 1022B 03 

4496. 

6 

0.6C22E 05 

0. 546000E 03 


0.9210E 01 

0.1114E 03 

5044. 

7 

0.6447E 05 

0. 47-0CC0E 03 


0.7724E 01 

0.1191E 03 

5514. 

8 

0.6896E 05 

0.S160C0E 03 


0.8458E 01 

0. 1275E 03 

603C. 

9 

0.1116E 06 

0.580100E 04 


0. 813 IE 02 

0.2089E C3 

11831. 

10 

0. 1158E 06 

0. 6660COE 03 


0.7969E 01 

0.2168E 03 

12497. 

11 

0.1192E 06 

0.553000E 03 


0.6604E 01 

0.2234E 03 

13CS0. 

12 

0.1233E 06 

0 . 668000 E 03 


0.7963E 01 

0.2314E 03 

13718. 

13 

0. 1261E 06 

0. 5100C0E 03 


0.60S9E 01 

0.23755 03 

14228. 

14 

0. 171 CE 06 

0.8592C0E 04 


0.88211 02 

0. 32571 03 

2282C. 

15 

0. 1748E 06 

0.824000E 03 


0.80C8E 01 

0.3337E 03 

23644. 

16 

0. 1780E 06 

0.703CC0E 03 


0.6824E 01 

0. 3405E 03 

24347. 

17 

0.1819E 06 

0.389000E 03 


0.783«l. 01 

0.3483E C3 

25236. 

18 

0..436*>E J6 

0. 1077CJ1 04 


0.8341a 01 

C.3 r - '»1E n 3 ■ 

26 313. 

19 

0.2252E 06 

0.988000E 04 


0.7596B 02 

0. 4J26Z 03 

36193. 

20 

0.2290E 06 

0. 1 13300 E 04 


0.7775E 01 

0.4404E 03 

37326. 

21 

0. 2297b 06 

0. 222000E 03 


0.1256E 01 

0.4417B 03 

37546. 

22 

0.14371 04 

0.870000E 02 


0. 261SE 01 

0.26165 01 

87. 

23 

0.4819E 05 

0. 371100E 04 


0.8563E 02 

0.8824E 02 

3798. 

24 

0. 4961E 05 

C. 14C0C0E 03 


0.261SE 01 

0.9086E 02 

3938. 

25 

0.9217E 05 

0 . 504500E 04 


0.7866E 02 

0.1695E 03 

8983. 

26 

0.9347E 05 

0. 1790C0E 03 


0. 2528E 01 

0.1720E 03 

9162. 

27 

0.9669E 05 

0.452CC0E 03 


0.5480B 01 

0. 1775E 03 

9614. 

28 

. 0.9780E C5 

0.1580001 03 


0. 1908E 01 

0.1794E 03 

9772. 

29 

0.9962B 05 

0. 262CCCE 03 


0.3160E 01 

0.1 826 E 03 

10034. 

30 

0. 1C35E 06 

0. 5620C0? 03 


0.6768S 01 

0. 1894E 03 

10596. 

31 

O'. 1 1 06 B 06 

0. 1 C8400E 04 


0. 1347E 02 

0.2028E C3 

1168C. 

32 

0. 1699B 06 

0. 1CCCC0E 05 


0. 11611 03 

0.3189E C3 

21680. 

33 

0. 1826E 06 

0.2839 CCE 04 


0.264SE 02 

0. 3454E 03 

24519. 

34 

0.1985E 06 

0.376100E 04 


0.2906E 02 

0.3744E 03 

2628C. 

35 

0. 2120E 06 

0. 344400E 04 


0.2650E 02 

0.4009E 03 

31724. 

36 

C.2260E 06 

0. 382900E 04 


0.2934E 02 

0. 43035 03 

35553. 

37 

0.2297B 06 

0. 111600E 04 

- 

0.7203E 01 

0.4375E 03 

36669. 



8IIE0TM 
. 80BBBB 
( 

TIHE 

(BBS) 

TIHEZ 

STSTES 

RELIAB 

PZ 

LEVEL OF ST SI EH 
1 2 
P?Z(K,1) PFZ(K,2) 

AVAILABILITY (BO. 

3 4 

PFZ(K,3) PFZ(K,4) 

OF PROP. 
5 

PPZ (K,5) 

TANKS) 

6 

PPZ (K,6) 

RELIABILITY OP 
THRUST SYSTEM 
REtZ(K) R22Z (K ) 

1 

480.0 

1 .ocoooo 

0.000000 

0.000000 

0.000032 

0.000004 

0.002958 

0. 997006 

1.000000 

1.00G000 

3 

504.0 

1,000000 

0.000000 

0. OCOOOO 

C. 000033 

C.C00004 

0.003106 

0.996856 

1.000000 

1 .oooooo 

3 

793.0 

1.CC0000 

o.cooooo 

a. oooooo 

0.000052 

0.000010 

0. 00487 3 

0.995065 

1.000000 

1 . oooooo 

• 

816.0 

1.000000 

0.000000 

0.000000 

0.000054 

0.000011 

0.005020 

0.994915 

1.000000 

1.000000 

5 

817.0 

1 . ocoooo 

0.000000 

o.cooooo 

C. 000054 

o.oocaii 

0.005026 

0.994909 

1.000000 

1.000000 

< 

1153.0 

1.CC0CC0 

o.ooocao 

o.ccoooo 

C. 000076 

0.000021 

0.00707 4 

0.992828 

1.000000 

1.000000 

7 

1176.0 

1.000000 

0.000000 

0. oooooo 

0.000078 

0.000022 

0.007220 

0.992680 

1 .ocoooo 

1.000000 

8 

1656.0 

1.000000 

o.ooococ 

o.caocoi 

0.000110 

0.000043 

0.010139 

0.989707 

1.000000 

1.000000 

9 

1680.0 

1.000000 

o.ooocoo 

C. 000001 

0.000111 

0.000045 

0.010285 

0. 989559 

1.000000 

1.000000 

10 

1681.0 

1.000000 

0.000000 

0.000001 

C. 000111 

0.000045 

0.010291 

0.989552 

1.000000 

1.000000 

11 

2016.0 

1.000000 

o.ooocoo 

O.CCOCOI 

C. 000134 

0.000064 

0.012319 

0.987482 

1.000000 

1.000000 

13 

2040.0 

1 . ocoooo 

0.000300 

0.000001 

0.000135 

0.000066 

0.012464 

0. 987334 

1.CC0000 

i.oocccc 

13 

3 176. 0 

1. 100000 

). 000001 

o. coocoi 

C. 00015? 

0.COOO84 

C. 014490 

0 .985262 

1.000000 

1 .000*30 

1* 

3»00.3 

0. 49*998 

C. 000301 

G.C09U )1 

l. 000160 

0.0000 J1 

0.014634 

0.985114 

1.000000 

i. ocoooo 

15 

2401.0 

0.999409 

0.000001 

0.000001 

C. 000160 

0.000091 

0.014640 

0.985109 

1.0C000C 

1 . cccccc 

16 

3712.0 

0. 999998 

0.000001 

0.000002 

0.000180 

0.00C115 

0. 016507 

0.983195 

1.000000 

1.000000 

17 

2736.0 

0.999998 

0.000001 

0.000002 

0.000182 

0.000118 

0.016652 

0.983047 

1.000000 

1.000000 

18 

3072.0 

0.999998 

0.000001 

0.000002 

0.000205 

0.000148 

0.018661 

0.980984 

1 . ocoooo 

1.000000 

19 

3096.0 

0.999998 

o.coocoi 

0.GC0002 

C. 000206 

0.000150 

0.018804 

0.980837 

i. oooooo 

1.000000 

30 

3097.0 

0.999999 

0.000001 

0.000032 

0.000207 

0.000150 

0.018810 

0.980831 

1 . ocoooo 

i .coooao 

31 

352B.0 

0.999997 

0.000001 

a . 000003 

0.000236 

0.C00195 

0.021377 

0.97819C 

i.ocococ 

1.000000 

33 

3888.0 

0.999622 

0. 000003 

0.CC00C4 

C. 000373 

C. 000236 

0.023509 

0.975877 

1.000000 

1.000000 

33 

4296.0 

0.999591 

0.000003 

0.000005 

0.000403 

0.000287 

0.025916 

0. 973388 

1 . cccccc 

i.caccoc 

3« 

4776.0 

0.999555 

0.000004 

0.000006 

C. 000438 

C.C0C354 

0.028734 

0.970466 

1.000000 

1.000000 

35 

5064.0 

0.999135 

0.000005 

C. 000007 

0.000459 

C. 000398 

0.030418 

0.968717 

1.000000 

1.000000 

36 

5065.0 

0.999135 

0.000005 

0.000007 

0.000459 

0.000398 

0.030423 

0.968712 

1.000000 

1.CCG0C0 

37 

5376.0 

0.999061 

O.OOOCQ5 

0.000008 

0.000482 

C..C0C448 

0.032235 

0.966826 

1.000000 

1.000000 

38 

5712.0 

0.998978 

0.000006 

0.000009 

C. 000507 

0.000505 

0.034185 

0.964793 

1.000000 

1.C0C0CC 

39 

6528.0 

0.998763 

0.000007 

0.000011 

0.000560 

0.000657 

0. 038890 

0.959873 

1 . ocoooo 

1 . cococo 

30 

6672.0 

0.998723 

0.000007 

0.000012 

C. 000579 

C. 000685 

0.039716 

0.959007 

1 . ocoooo 

1.000000 

31 

6816.0 

0. 998682 

0.C0CCC8 

C.C00012 

0.000590 

0.000715 

0.040540 

0.958142 

1.000000 

1.0000CC 

33 

7200.0 

0.955839 

.0. 000008 

0.000014 

0.000619 

0.000796 

0.042732 

0.955839 

1.000000 

1.000000 

33 

7536.0 

0.953826 

0.000CC9 

0.000015 

C. 000645 

0.000871 

C. 044642 

0.953826 

1.000000 

1.000000 

39 

7537.0 

0.953821 

0.C00CC9 

0.000015 

0.000645 

0.000871 

0.044647 

0.953821 

1 . ocoooo 

1.000000 

35 

7872.0 

0.951819 

0.000010 

0.000016 

0.000671 

0. 000948 

0.046544 

0.951819 

0.999999 

0.999999 

36- 

7873.0 

0. 951814 

0.000010 

C.C00C16 

C. 000671 

0.000949 

0.046549 

0.951814 

0.999999 

0.999999 

37 

8208.0 

0.949815 

0.000011 

0.000018 

0.000697 

0.001029 

0.048439 

0.949815 

0.999999 

0.999999 
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MANEUVER 

MUHBEB 

K 

MUHBEB OP 
CICI.ES 

CICL1Z CTCL2Z 

PILL 

VALVES 

BPVZ 

HI PBESS 
TASK 
BTZ 

PROP TASK 
HELD SIMS 
8PTZ 

PBOP 

VALVE 

BSP1Z 

TAMK 

SYSTEH 

BSP2Z 

BLADDER 
BELIASILITT 
RBI Z RB2Z 

PBESS OPE 
TRANSDUCES 
BPD1Z BPD2Z 

1 

1 

1 

0 

1.000000 

0.999999 

0.999999 

1.000000 

1.000000 

0.999520 

0.999520 

0.999998 

1.000000 

2 

2 

0 

1.000000 

0.999999 

0.999999 

1. 000000 

1.000000 

0.999496 

0.999496 

0.999995 

1.000000 

3 

3 

0 

1.000000 

0.999998 

0.999998 

1.000000 

1.000000 

0.999208 

0.999208 

0.999992 

1.000000 

• 

4 

0 

1.000000 

0.999998 

0.999998 

1.000000 

1.000000 

0.999184 

0.999184 

0.999960 

1.000000 

5 

5 

0 

1.000000 

0.999998 

0. 999998 

1.000000 

1.000000 

0.999183 

0.999183 

0.999956 

1.000000 

6 

6 

C 

1.000000 

0.999997 

0.999997 

1.000000 

1 .000000 

0.998849 

0.998849 

0.999951 

1.000030 

7 

7 

C 

1.000000 

0.999996 

0.999996 

1.000000 

1.000000 

0. 9^8825 

0. 999825 

C. 999946 

1 . ccoooo 

8 

8 

0 

1.000000 

0. 999995 

0. 999995 

1.COCOCC 

1.000000 

0.388335 

0.499345 

J. 999941 

1.000003 

9 

9 

c 

1.0C0000 

0.999995 

0.999995 

1.00 30 00 

1 .OJOOOj 

0. ?98 J21 

0. 398321 

0. 999395 

1. 000330 

10 

10 

c 

1. 000300 

0.999995 

0.999995 

1.000000 

1.000000 

C. 098320 

0.999320 

0.999879 

1.000030 

11 

11 

0 

0.999999 

0.999994 

0.999994 

1 . cococo 

1.000000 

0.997986 

0.997986 

0.999373 

1.000000 

12 

12 

0 

0.999999 

0.999994 

0.999994 

1.000000 

1.000000 

0.997962 

0.997962 

0.999867 

1.000300 

13 

13 

0 

0.999999 

0.999993 

0.999993 

1.000000 

1.000000 

0.997627 

0.997627 

0.999862 

1.000000 

14 

14 

0 

0. 999999 

0.999993 

0.999993 

1.000000 

1 .000000 

0.997603 

0.997603 

0. 999778 

1.000000 

15 

15 

0 

0.999999 

0.999993 

0.999993 

1.000000 

1.000000 

0.997602 

0.997602 

0.999770 

1 .ocoooo 

16 

16 

0 

0.999999 

0.999992 

0.999992 

1.000000 

1.000000 

0.997292 

0.997292 

0. 999763 

1.000000 

17 

17 

0 

0.999999 

0.999992 

0.999992 

1.000000 

1 .000000 

0.997268 

0.997268 

0.999755 

1.000000 

18 

18 

0 

0.999998 

0.999991 

0.999991 

1.000000 

1.000000 

0.996933 

0.996933 

0.999744 

1 . ocoooo 

19 

19 

0 

0.999998 

0.999991 

0. 999991 

1.000000 

1.000000 

0.996909 

0.996909 

0.999648 

1.000000 

20 

20 

0 

0.999998 

0.999991 

0.999991 

1.000000 

1.000000 

0.996908 

0.996908 

0.999637 

1.000000 

21 

21 

0 

0.999998 

0.999989 

0.999989 

0. 999999 

0.999999 

0.996478 

0.996478 

0.999635 

1.000000 

22 

21 

1 

0.999997 

0.999988 

0.999988 

0.999999 

0.999999 

0. 996119 

0.996120 

0.999635 

0.999999 

23 

21 

2 

0.999997 

0.999987 

0.999987 

0.999999 

0.999999 

0.995713 

0.995713 

0.999635 

0.999963 

24 

21 

3 

0.999996 

0.999986 

0.999986 

0. 999999 

0.999999 

0.995235 

0.995235 

0.999635 

0.999962 

25 

21 

4 

0.999995 

0.999985 

0.999985 

0.999999 

0.999999 

0.994948 

0.994949 

0.999635 

0.999913 

26 

21 

S 

0.999995 

0.999985 

0.999985 

0.999999 

0.999999 

0. 994947 

0. 994546 

0.999635 

0.999911 

27 

* 21 

6 

0.999995 

0.999984 

0.999984 

0.S99999 

0.999999 

0.994638 

0.994638 

0.999635 

0.999907 

28 

21 

7 

0. 999994 

0.999983 “ 

0.999983 

0.999998 

0.999998 

0.994304 

0.994304 

0.999635 

0.999905 

29 

21 

8 

0.999992 

0.999980 

0.999980 

0.999998 

0.999998 

0.993493 

0. 993493 

0.999635 

C. 999902 

30 

21 

9 

0.999992 

0.999980 

0.999980 

0. 999998 

0.999998 

0.993350 

0.993350 

0.999635 

0.999897 

31 

21 

10 

0.999992 

0.999.980 

0.999980 

0.999998 

0.999998 

0.993207 

0.993207 

0. 999635 

0.999886 

32 

21 

12 

0.999991 

0.999978 

0.999978 

0. 999997 

0.999997 

0. 992826 

0. 992826 

C. 995635 

0.9S9789 

33 

21 

13 

0.999990 

0.999977 

0.999977 

0. 999997 

0.999997 

0.992492 

0.992492 

0.999635 

0.999762 

34 

21 

14 

0.999990 

0.999977 

0.999977 

0.999997 

0.999997 

0.992491 

0. 992491 

0.599635 

0.999725 

35 

21 

15 

0.999989 

0.999976 

0.999976 

0. 999997 

0.999997 

0.992159 

0.992159 

0.999635 

0. 999692 

36 

21 

1 € 

0.999989 

0.999976 

0.999976 

0. 999997 

C. 999997 

0.992158 

0.992158 

0.999635 

0.999654 

37 

21 

17 

0.999988 

0.999975 

0.999975 

0.999997 

0.999997 

0. 991825 

0.991826 

0.999635 

0.999644 
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» 

HABE0YI8 FILTER THBDSTER THRUSTER HALF-SYS HELD CONNECTION SELIA3 

NOBEEH ASSEMBLY VALYE SYSTEM RELIABILITY CON. V&LV HP SYS PRO SYS THR SYS 

K RF1Z RF2Z SST1Z RST2Z RTH1Z BTH2Z BZ RLHPZ RLPPZ RLTZ 

1 0.999999 0.999999 1.000000 1.000000 0.999991 0.999991 1.000000 0.999986 0.999986 0.999986 

2 0.999999 0.999999 1.000000 1.000000 0.999991 0.999991 1.000000 0.999985 0.999985 0.999965 

3 0.999998 0.999998 1.000000 1.000000 0.999986 0.999986 1.000000 0.999976 0.999976 0.999976 

4 0.999998 0.999998 1.C00000 1.000000 0.999985 0.999985 1.000000 0.999976 0.999976 0.999976 

5 0.999998 0.999998 1.000000 1.000000 0.999985 0.999985 1.000000 0.999976 0.999976 0.999976 

6 0.999997 0.999997 1.000000 1.000000 0.999979 0.999979 1.000000 0.999965 0.999965 0.999965 

7 0.999996 0.999996 1.000000 1.000000 0.999979 0.999979 1.000000 0.999965 0.999965 0.999965 

8 0.999996 0.999995 1.000000 1.000000 0.999970 0.999970 1.000000 0.999950 C. 999950 0.999950 

9 0.999996 0.999995 1.000000 1.000000 0.999969 0.999970 1.000000 0.999950 0.999950 0.999950 

10 0.999996 0.999995 1.000000 1.000000 0.999969 0.999970 1.000000 0.999950 0.999950 0.999950 

11 0.999995 0.999994 1.000000 1.000000 0.999963 0.999964 1.000000 0.999940 0.999940 0.999940 

12 0.999995 0.999994 1.000000 I.COOOCO 0.999963 0.999963 1.000000 0.999939 0.999939 0.999939 

13 0.999994 0.999993 1.000000 1.000000 0.999957 0.999957 1.000000 0.999929 0.9S9S29 0.999929 

14 0.999995 0.999993 1.000000 1.000000 0.999956 0.999957 1.000000 0.999928 0.999928 0.999928 

15 0.999995 0.999993 1.000000 1 . 000000 0.999956 0.999957 1.000000 0.999928 0.999928 0.999928 

16 0.999994 0.999992 1.000000 1.000000 0.999950 0.999951 I.OCOOOO 0.999919 0.999919 0.999919 

17 0.999994 0.999992 1.000000 1.000000 0.999950 0.999951 I.OCOOOO 0.999918 0.999918 0.999918 

18 0.999993 0.999991 1.00000C 1.000000 0.999944 0.999945 1.000000 0.999908 0.999908 0.999908 

19 0.999994 0.999991 1.000000 1.000000 0.999943 0.999944 1.000000 0.9999C7 0.999907 0.999907 

20 0.999994 0.999991 1.000000 I.OCOOOO 0.999943 0.999944 1.000000 0.999907 0.999907 0.999907 

21 0.999992 0.999989 0.999999 0.999999 0.999935 0.999937 0.999999 0.999894 C. 999894 0.999894 

22 0.999991 0.999988 0.999999 0.999999 0.999929 0.999930 0.999999 0.999883 0.999883 0.999883 

23 0.999990 0.999987 0.999999 0.999999 0.999921 0.999923 0.999999 0.999871 0.999871 0.999871 

24 0.999988 0.999985 0.999999 0.999999 0.999913 0.999914 0.999999 0.999857 0.999657 C. 999857 

25 0.999988 0.999984 0.999999 0.999999 0.999907 0.999909 0.999999 0.999848 0.999848 0.999848 

26 0.999988 0.999984 0.999999 0.999999 0.999907 0.999909 0.999999 0.999848 0.999848 0.999848 

27 0.999987 0.999983 0.999999 0.999999 0.999902 0.999903 0.999999 0.999839 0.999839 C. 999839 

28 0.999986 0.999982 0.999998 0.999998 0.999896 0.999897 0.999998 0.999829 0.999829 0.999829 

29 0.999983 0.999980 0.999998 0.999998 0.999881 0.999882 0.999998 0.999804 0.999804 0.9998C4 

30 0.999983 0.999979 0.999998 0.999998 0.999878 0.999879 0.999998 0.999800 C. 999800 C. 999800 

31 0.999982 0.999979 0.999998 0.999998 0.999876 0.999877 0.999998 0.999796 0.999796 0.999796 

32 0.999981 0.999977 0.999997 0.999997 0.999869 0.999870 0.999997 0.999784 0.999784 0.999784 

33 0.999980 0.999975 0.999997 0.999997 0.999863 0.999863 0.999997 0.999774 C. 999774 0.999774 

34 0.999980 0.999975 0.999997 0.999997 0.999863 C. 999863 0.999997 0.999774 0.999774 0.999774 

35 0.999979 0.999974 0.999997 0.999997 0.999857 0.999857 0.999997 0.999764 0.999764 C. 999764 

36 0.999979 0.999973 0.999997 0.999997 0.999857 C. 999857 0.999997 0.999764 0.999764 0.999764 

37 0.999978 0.999972 0.999997 0.999997 0.999851 0.999851 0.999997 0.999754 0.999754 0.999754 
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APPENDIX C: Valve Configurations . ■ 

- One of the weakest links in auxiliary propulsion systems is the 
solenoid valve. Four commonly used valve redundancy configurations are: dual 
series, dual parallel, quad, and quad connected. These are illustrated below 
along with equations for their reliability. Two basic failure modes are con- 
sidered; i.e., open failures and closed failures. Both an open and a closed 
failure will lead to a system failure in the absence of redundancy. 

In general, let 

R q = probability of no open failure 

R c = probability of no closed failure 

Q q = 1-R q = probability of open failure 

Q c = 1-R^ ~ probability of closed failure 

Q g = Q q + Q c = probability of failure of a single valve 

R = valve system reliability 

Q = valve system unreliability 



R . 1 . q s . 1 - q o . q c . 1 - (1 - R„) - U - R c ) 

Dual Series 

The various possible outcomes are shown below 


t = R + R - 1 
° c 




Therefore, 



The various possible outcomes are shown below 



Valve # I 
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A system failure will occur if either valve fails in the open position or if 
both valves fail in the closed position. Therefore, 


R = 1 - 1 - 


2 2 
R = R - (1 - R ) 
o c 


a - Q o ) 2 + Q c 2 




2 2 
where 1 - (1 - Q q ) is the probability of an open failure and Q c is the 


probability of a closed failure. 

r-C&J— i 


Quad 




The various possible outcomes can be structured in a manner similar to 
that shown above. It should be noted that the probability of a system failure is 
the sum of the probability of a failure due to an open valve arrangement and that 

I 

due to a closed valve arrangement. Therefore, 


probability of closed failure 


- [i - a - q c > 2 _ 


2,2 


probability of open failure = 1 - (1 - Q ) 


R 


■ i - f 1 - <i - v 2 ] - 

= [ 1 ; a- R /] 2 - [i-k c 2 ] 


1 + (1 - Q q 2 ) 2 


r-t^3-pC^0-i 

Quad Connected 

H&H-c&h-J 

In a manner similar to the Quad above, 
probability of closed failure = 1 - £l - Q c 2 J 



/ 
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probability of open failure = jJL — (1 — Q q ) j ^ 

R - 1 - [l - <1 -Q 0 ) 2 ] 2 - 1+ [l -Q c 2 ] 2 
R- [l- U-R c ) 2 ] 2 - [l-R„ 2 ] 2 



